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Abstract 

A generalization of the Gell-Mann-Low Theorem is applied to bound state cal- 
culations in Yukawa theory. The resulting effective Schrodinger equation is solved 
numerically for two-fermion bound states with the exchange of a massless boson. The 
complete low-lying bound state spectrum is obtained for different ratios of the con- 
stituent masses. No abnormal solutions are found. We show the consistency of the 
non-relativistic and one-body limits and discuss the special cases of identical fermions 
and fermion-antifermion states. To our knowledge, this is the first consistent calcula- 
tion of bound states in pure Yukawa theory (without UV cutoff). 



1 Introduction 

During the more than 50 years that have passed since Bethe and Salpeter formulated their 
famous equation [UI2], the calculation of relativistic bound states has proved to be one of 
the truely hard problems in quantum field theory. For most of the numerous proposals for 
its solution, a model theory of two scalar bosons interacting via the exchange of a third 
scalar has served as a first testing ground. In the case of the Bethe-Salpeter equation itself, 
it was shown by Wick and Cutkosky that this model theory has an analytical solution in 
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the popular "ladder approximation" to the equation, in case that the exchanged boson is 
massless |3;. It is hence quite surprising that for one of the most natural generalizations of 
this model where two spin 1/2 fermions interact through the exchange of a scalar boson, not 
a single consistent formalism has been devised for its solution to date. 

In this introduction, we will give a brief review of the most important intents to calculate 
relativistic two-fermion bound states in this model, i.e., in Yukawa theory. It was clear even 
before Wick and Cutkosky's solution of the purely scalar model that the case of fermionic 
constituents means a lot more to the Bethe-Salpeter equation than just a technical compli- 
cation due to the inclusion of spin degrees of freedom jlj: the fact that the kernel of the 
equation in the ladder approximation (after a Wick rotation) is not of Fredholm type in this 
case represents a fundamental difficulty for analytic as much as for numerical investigations. 
In the case of equal constituent masses and zero boson mass, it was shown by Goldstein |H] 
(later corrected and improved upon by Green jH]), originally for the subspace of pseudoscalar 
bound states only, that the ladder approximation gives a continuum of coupling constants 
corresponding to a "tightly bound" state with zero mass, i.e., such that the binding energy 
completely compensates the masses of the constituents, while one expects a series of discrete 
values if a discrete spectrum of energy eigenvalues is to exist. The definite version of the 
argument is reviewed in Ref . [7] (see the references there) . In the case of Yukawa theory, 
this is essentially all that is known about the solutions of the Bethe-Salpeter equation. 

The exchange of a scalar boson between spin 1/2 fermions has come to play an impor- 
tant role in the description of the nucleon-nucleon interaction in the context of one-boson 
exchange (OBE) models. Although the scalar boson exchange is used as an effective and 
approximate description of an actual two-pion exchange in this case, it is believed to give the 
dominant attractive contribution to the intermediate-range potential. In the numerous nu- 
merical calculations for the deuteron system within OBE models, the scalar boson exchange 
is always accompanied by the exchange of other mesons, most importantly pions and rho and 
omega mesons. For this type of calculations, the Bethe-Salpeter equation as well as several 
of its three-dimensional reductions have been employed, notably the Blankenbecler-Sugar- 
Logunov-Tavkhelidze (BSLT) equation [S] and the Gross or spectator equation 0. The 
latter equations avoid problems with abnormal solutions, the one-body limit, and others 
that arise in the Bethe-Salpeter equation (in the ladder approximation) jHEII], see, however, 
Ref. [TT] . Among the vast literature on OBE models, we mention the work of Tjon and 
collaborators ^2] who use the Bethe-Salpeter equation, the work of Gross and collaborators 
[T3] who base their analysis on the Gross equation, and the quite complete "Bonn model" by 
Holinde, Machleidt, et al. reviewed in Ref. |14| . which employs the Bloch-Horowitz scheme 
[To] for the bound state calculations, a formalism that is not based on the Bethe-Salpeter 
equation. 

In all these approaches, a momentum cutoff is introduced in order to make the equations 
well-defined. It is interpreted phenomenologically as representing the spatial extension of 
the nucleons and mesons. Although we have not found a corresponding calculation in the 
published literature, it appears plausible in view of the singular nature of the Bethe-Salpeter 
kernel in the ladder approximation that the introduction of a momentum cutoff would also 
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be necessary in the case of a pure Yukawa interaction, and that the corresponding numerical 
results depend on the value of the cutoff. Note, however, the work of Gari and collabo- 
rators ^H] who dispense with a momentum cutoff and rather determine the nucleon and 
meson structures self-consistently through loop corrections within the same nucleon-meson 
model. It is perhaps not a coincidence that their approach employs the Okubo transforma- 
tion technique ^7j, which is unrelated to the Bethe-Salpeter equation but very similar to 
the formalism proposed in the present paper. 

More recently, calculations of bound states in Yukawa theory have been attempted in 
light-front quantization, both in a Tamm-Dancoff approximation jTH] and in the so-called 
covariant light-front dynamics |X§]. Just as in the equal-time quantized theory, an additional 
momentum cutoff has to be introduced for the one-boson exchange, and the solutions for the 
bound states turn out to depend on this cutoff. It was pointed out by van Iersel and Bakker 
[2*U] that the light front ladder approximation is incomplete as long as instantaneous terms 
are not taken into account. One possibility to include such terms was presented in Ref. |2I] 
and may help to solve the renormalization problem in the future. 

In this paper, we propose to use a generalization of the Gell- Mann-Low theorem [2] 
for the calculation of relativistic bound states. After the obligatory application to the purely 
scalar model in a previous publication [23], we now turn to Yukawa theory. Surprisingly, 
the application of the generalized Gell-Mann-Low theorem turns out to be straightforward 
and, except for the necessary technical complications due to the spin degrees of freedom, 
completely analogous to the scalar case. No inconsistencies whatsoever have been found, 
and no necessity for the introduction of a cutoff (for the interaction between the fermions) 
arises. For the numerical calculation of the bound state spectrum, we will focus on the case 
of a massless exchanged boson, for the following reasons: first, the massless case is the most 
singular one, and if the formalism works in this case, it is expected to be applicable to the 
case of a massive exchanged boson as well. Second, the solutions in the non-relativistic limit 
are known analytically in this case (Coulomb potential) and our numerical solutions can be 
compared against them. In particular, one would like to find the degeneracies characteristic 
of the non-relativistic limit. Finally, we plan to extract the lowest-order fine and hyperfine 
structure from our effective Hamiltonian in the near future, and this can be done analytically 
only in the massless (Coulomb) case. 

The organization of the paper is as follows: in the next section, we introduce the gen- 
eralization of the Gell-Mann-Low theorem on which our approach is based, and comment 
on aspects relevant to the present work. We also present, in the same section, the effective 
Hamiltonians in the zero-, one-, and two-fermion sectors which result from the application 
of the generalized Gell-Mann-Low theorem to lowest non-trivial order. They describe the 
renormalization of the vacuum energy, the fermion mass renormalization, and the effective 
potential for the two-fermion dynamics, respectively. In Section 3, we perform a series of 
non-trivial analytical checks on the results obtained in Section 2. In particular, we investi- 
gate the non-relativistic and one-body limits of the effective Schrodinger equation, replace 
fermions by antifermions, and consider the case of identical constituents. All the properties 
of our effective Hamiltonian description turn out to be in accord with physical expectations. 
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We present numerical bound state solutions in Section 4, for fine structure constants be- 
tween zero and one and arbitrary ratios of the constituent masses. The eight lowest-lying 
states are calculated in each case, including states with non-zero relative orbital angular 
momentum and mixtures of spin singlet and triplet states. The non-relativistic, one-body, 
and equal-mass limits are discussed in detail. Section 5 contains our conclusions. There are 
also four appendices. Appendices A and B are concerned with the appearing loop corrections 
and their regularization in different schemes from a Hamiltonian, not manifestly covariant 
perspective. They prepare the ground for future calculations at higher orders in the expan- 
sion of the effective Hamiltonian, where regularization and renormalization are expected to 
become central issues. Appendices C and D present several formulas which are used in the 
analytical separation of angular and spin variables in the effective Schrodinger equation. 

2 The Bloch- Wilson Hamiltonian 

The physical idea behind the Bloch- Wilson or effective Hamiltonian is very similar to the 
Born-Oppenheimer approximation: the integration over the "light" degrees of freedom gen- 
erates the dynamics for the "heavy" degrees of freedom (the constituents in our case), where 
the "light" degrees of freedom are given by the so-called interacting "virtual clouds" around 
the constituents. The virtual clouds are, in turn, created by the constituents as described 
here through an adiabatic process. Technically, a generalization of the Gell-Mann-Low the- 
orem is used which we will briefly describe in the following (for further details, see Ref. 



The adiabatic evolution operator U e provides a map from an eigenstate of the free Hamil- 
tonian to an exact eigenstate of the full interaction Hamiltonian, if its application to the free 
eigenstate, after a suitable normalization, is well defined. This is the content of the original 
theorem by Gell-Mann and Low 0. To this end, the full Hamiltonian H is split into a free 
part H and an interacting part H\. The adiabatic damping is implemented by replacing 



An eigenstate \<j>) of H = lim^-oo H (t) is then evolved according to the time-dependent 
Schrodinger equation with Hamiltonian H(t) from t — ► — oo to the state U e \4>) at t = 0. 
According to the adiabatic theorem, in the limit e — > the state U e \4>) is an eigenstate of 
H(t = 0) = H. A possible infinite phase (see, for example, [21]) is taken care of by imposing 
the normalization condition 



22). 



H = H + H\ through 



H{t) = H + e 



t\ H , 




(2) 



so that 




(3) 
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As a generalization of this theorem, it can be shown that the adiabatic evolution 
operator also maps i7 -invariant subspaces fl , H Q Q C fi , to //-invariant subspaces fl, 
HQ C Q, if its application to the HQ-invariant subspace, after a suitable normalization, is 
well defined. When a subspace is invariant under a hermitian operator, this operator can 
be diagonalized in the subspace, in other words, the subspace is a direct sum of eigenspaces 
of the operator, in this case the Hamiltonian. Since the only important property of the 
adiabatic evolution operator in this context is its iZ-invariant image, there is still an infinity 
of possible "normalizations" of the operator that map between the same subspaces Qq and 
Q. Here we choose the normalization condition 

-Po^BW — Ifio > (4) 

where P is the orthogonal projector to Q and 1q the identity operator in Q . Ubw is the 
"Bloch- Wilson operator", the normalized version of the adiabatic evolution operator given 
explicitly by 

Ubw = UMUJ- 1 . (5) 

Eq. naturally generalizes the normalization condition (j2J) in the original Gell-Mann-Low 
theorem. The subspace Qq has to be small enough in order that it is possible in practice to 
diagonalize the Hamiltonian in Q = Ubw^o, at least numerically. On the other hand, it has 
to be large enough for the Bloch- Wilson operator to be well defined. This is in principle a 
subtle issue and may depend on the normalization condition chosen. We will limit ourselves 
here to show that everything is well defined for a natural choice of the subspace (for the field 
theory and the perturbative order considered here, see the discussion below). In the case of 
a two- constituent bound state, a natural choice for the if -i nvar i an t subspace Qq is the space 
of all states of the two constituents as free particles. This will be mapped by the normalized 
adiabatic evolution operator to the .//-invariant subspace Q, a direct sum of eigenstates of 
the full Hamiltonian, which is expected to coincide with the space of all physical two-particle 
states, scattering as well as bound states. 

A further simplification is obtained by realizing that the Bloch- Wilson operator effects 
a similarity transformation between the subspaces Qq and Q. Since the free subspace Qq is 
usually much easier to work in, it is convenient to similarity transform the full Hamiltonian 
back to this subspace, thereby defining an effective or Bloch- Wilson Hamiltonian Hbw i n 
Q , 

Hbw = PqHUbw , (6) 

where we have taken into account that Pq\q is the inverse map to Ubw as a consequence of the 
normalization (j3J). In the above example, the effective Hamiltonian acts on the subspace of all 
states of the free constituents. Its eigenvalues are exactly equal to the eigenvalues of the full 
Hamiltonian (in the corresponding subspace Q), and there is a one-to-one relation between 
the corresponding eigenstates. Mathematically, the time-independent Schrodinger equation 
for the effective Hamiltonian takes the form of a non-relativistic Schrodinger equation for the 
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two constituents, with the relativistic expression for the kinetic energies and a complicated 
(generally non-local and non-hermitian) interaction term. With the normalization chosen, 
the eigenstates \<f>) of the effective Hamiltonian have the immediate meaning of wavefunctions 
for the constituents: 

|0>=Po(%w|0», (7) 

hence \<p) is the two-particle component of the complicated exact eigenstate of H, Ubw\<P)- 
Obviously, the Bloch- Wilson operator Ubw cannot be determined exactly in any in- 
teresting practical application. However, the adiabatic evolution operator U e , and hence 
Ubw — U e (PoU e )~ l , have a well-known perturbative expansion, the Dyson series: 

00 (—j\ n f° f° 
Ue = Y,-^~ / / dt n e-<^ + - + ^T[H 1 (t 1 )---H 1 (t n )}, (8) 

where 

H 1 (t) = e iHot H 1 e- iHot . (9) 

Every order in this expansion determines an approximation to the effective Hamiltonian 
and consequently of its eigenvalues and eigenstates. To compare with, in the Bethe-Salpeter 
approach an infinite part of the Dyson series has to be summed up to find an approximation to 
the bound states. The sum is performed in practice implicitly by solving an integral equation 
(the homogeneous Bethe-Salpeter equation) with a perturbative approximation to its kernel, 
and the counterpart in our Hamiltonian approach is the solution of the (approximated) 
effective Schrodinger equation. All of the above is exemplified in the application to a simple 
scalar model in j2S]. 

Finally, let us comment on the sense in which the adiabatic evolution operator, or rather 
its normalized version Ubw-, i s well-defined in Q . In the proof of the generalized Gell-Mann- 
Low theorem, as well as in the proof of the original theorem, the adiabatic evolution operator 
is treated as a formal power series (in Hi) given by the Dyson series. Consequently, Ubw 
is considered well-defined if every term of the corresponding series is well-defined (finite) 
disregarding convergence properties of the series as a whole. However, the UV divergencies 
of the usual covariant perturbation theory cannot be escaped: they appear in the context of 
the Bloch- Wilson Hamiltonian (in the most straightforward approach) as divergencies of the 
three-dimensional momentum integrals for large momenta. In the scalar model considered 
before it was shown that these divergencies can be dealt with (to the lowest non-trivial order 
considered there) rather easily, and, in fact, in the same way as in covariant perturbation 
theory. On the other hand, divergencies of the IR type are considered to potentially invalidate 
the existence of the Bloch- Wilson operator (for the specific subspace Qq chosen) . In the time- 
independent version presented in [221, which is obtained by performing the time integrals in 
the Dyson series, the IR divergencies arise from vanishing energy denominators. They can 
usually be avoided by appropriately enlarging the subspace Qo. However, the usefulness of 
the entire approach rests on the possibility of choosing a subspace f2 which is manageably 
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small, ideally the space of all states of the two constituents as free particles in the example 
mentioned above. Here, as for the scalar model before, we take a pragmatic approach: we 
will be content with the fact that no IR divergence arises to the perturbative order and for 
the choice of Qq considered. 

We will now apply these ideas to Yukawa theory. To simplify matters, we consider 
two different fermion species, A and B, which interact with a scalar boson. Hence the 
Hamiltonian of the theory consists of the free Hamiltonians for the fermions A and B and 
the scalar boson, and the (normal-ordered) interaction term 



H 1 = g J d 3 x: [Vu(x)^(x) + ip B (x)if) B (x)] <p(x) : 



(10) 



We are interested in bound states of one fermion A and one fermion B. Other bound states, 
of identical fermions or fermion and antifermion (e.g., in a theory with only one fermion 
species), are simply related to bound states of type AB and can be treated by exactly the 
same method. They will be discussed in some detail in the next section. The application 
of the generalized Gell-Mann-Low theorem to lowest non-trivial order follows in very close 
analogy the purely scalar case [2H1- In particular, for a subspace fl with a fixed number 
of (free) particles PqHiPq = 0, and the effective Hamiltonian becomes to lowest non-trivial 
order 

H BW = H P - i I dte^PoH^H^Po + 0(H*) (11) 

J — oo 

(the limit e — > being understood). For the rest of this chapter, we will discuss in detail 
the cases where Qq is the subspace of zero, one (fermionic), and two (AB) particles. We 
will go through the essential steps briefly, present the corresponding results, and take the 
opportunity to comment on several issues that have not received due attention in [23]. For 
illustrative purposes, it may be helpful to have a look at the diagrams presented there which 
are the same in the present case, only that spin labels have to be added to the external 
fermion lines. 

2.1 Vacuum state and zero-point energy renormalization 

Beginning with the vacuum state, we consider the Bloch- Wilson Hamiltonian for the subspace 
Q = C|0), or, equivalently, for the orthogonal projector on Q , P = |0)(0|. In this case, the 
application of the generalized version of the Gell-Mann-Low theorem is equivalent to the 
original Gell-Mann-Low theorem. The result for the Bloch- Wilson Hamiltonian is -ffew — 
EyPo, where (by Wick's theorem) 

Ey — Eq 

= ig 2 j dt I d 3 x d 3 x'A F {0 - t, x - x') tr [£#(0 - t, x - x.')S$(t - 0, x' - x)] 

J — oo J 

+ ig 2 j dte~ e W f d 3 xd 3 x'A F (0-t,x-x')tr [£f (0 - t, x - x')S#(t - 0, x' - x)] . (12) 

J — oo J 
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Here Eq is the vacuum energy of the free theory and Sp B and are the fermionic (for 
particles A, B) and the bosonic Feynman propagators. In our conventions, 



#(«,«) = « /t& r + Z A «-*■-"' (13) 

(27r) 4 — m A + «e 



d 3 p 



(2n) 3 2E£ L 



0(t) (p • 7 + m A ) e" i(B p '~ p ' x) - 9{-t) (p • 7 - m A ) e i(s p'- p - x) 



(14) 



in the covariant and non-covariant representation, respectively. We are using the following 
shorthands for the relativistic kinetic energies: 



E* = ^m\ B + p 2 , o; p = Vi" 2 + P 2 (15) 

(fi is the boson mass). The bosonic propagator Ap(t, x) is simply given by the m^-coefficient 
of Sp(t,x.) (and the substitution tjia — > //). In Eq. (|T2|) . note the opposite sign as compared 
to the bosonic case [23]. 

We can use time translation invariance of the integrands on the r.h.s. of Eq. ()12j) to 
move the vertices from times t < and t = to t = and — t > 0. Mathematically, this 
corresponds to replacing t by — t and exchanging x and x', just as we did for the one-particle 
states in the scalar case [23]. If we, furthermore, use the invariance of the bosonic propagator 
under the change of the direction of propagation and apply these manipulations to half the 
r.h.s. of Eq. (|12jl . we arrive (in the limit e — > 0) at the usual expression for the one-loop 
correction to the vacuum energy in covariant perturbation theory [23] . 
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J d A x d A x' 5(x°)A F (x - x') tr [S${x - x')Sj{x' - 
+ f d4 xd 4 x'5(x°)A F (x-x')tr [Sf(x - x')S§(x' - x)] . (16) 
The appearance of S(x°) eliminates a factor 

POO 

2tt5(0) = dt (17) 



00 



(due to time translation invariance) in the result, and leaves us with 

(2tt) 3 (5 (3) (0) = / d 3 x (18) 



(due to spatial translation invariance), to be interpreted as the volume of space. 

We can immediately express Eqs. (fT2*j) and ()16j) in momentum space, by use of the non- 
covariant and covariant expressions for the Feynman propagators in momentum space, Eqs. 
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(JHJ) and (JT3J), respectively. The results are 



Ey — Eq 



d 3 p d 3 p' 



1 4(-£ p ^ + p-p' + m 2 4 ) 



+ 9 2 



2 



g 2 
+ — 

2 



(2tt) 3 (2tt) 3 2£^ 2E$ 2u p+p , E£ + E$ + u p+p , 
d 3 p d 3 p' 1 4 (-E%E$ + p • p' + m 2 B ) 



(2tt) 3 (2tt) 3 2EZ IE* 2u 
d 4 p d 4 p' 



p+p' 



4(p • p' + m 2 ^) 



p+p' 



(2tt) 3 <5 (3) (0) 



m\ + ie] 



(2vr) 4 (2vr) 4 [(p - p') 2 - A* 2 + *e][p 2 - + ie][p /2 
(i 4 p d 4 p' 4(p ■ p' + m B ) 

(2tt) 4 (2vr) 4 [(p - pi) 2 - /i 2 + ie] [p 2 - m\ + ie] [p' 2 - m? B + ie] 



(2tt) 3 (5 (3) (0). 



(19) 



The equivalence between the two expressions can, of course, be established directly by inte- 
grating over p and p' in the manifestly covariant expression. However, there is a subtlety 
involved in the integration which may be of interest for future calculations at higher orders 
and will be described in detail in Appendix [XJ The results (fTTJj) are highly UV divergent 
and will be treated as formal expressions only. It is, however, important to remark that 
they are IR finite in the limit fi — > (we are here only interested in massive constituents 
m A ,m B ^ 0). 



2.2 One-fermion states and mass renormalization 

We will now turn to the one-fermion states, considering for concreteness fermions of type A. 
The if -invariant subspace considered is, correspondingly, 

Vl = span { \p A , r) \p A e M 3 , r = 1, 2} , (20) 

where \p A ,r) stands for a state of one fermion A with 3- momentum p^ and spin orientation r 
(we do not fix the basis to be used in spin space yet). The one-fermion states are normalized 
in a non-covariant fashion, 

(Pa, r\p' A , s) = (2n) 3 5(p A - P A ) $rs ■ (21) 
The application of the generalized Gell-Mann-Low theorem to lowest non-trivial order 
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yields the following matrix elements of the Bloch- Wilson Hamiltonian in Q , 
(p A ,r\H BW \p' A ,s) = {E V + E^ A ) (2ir) 3 6(p A -p' A )6 rs 



dte~ eltl / d 3 xd 3 x' 



dte~ eltl / d 3 xd 3 x' 



^(0, x)^(0 - t, x - x')^,Jt, x') A F (0 - t, x - x') 



$£ Atr (t,x!)S$(t - 0,x' - x)^ )S (0,x) A F {t - 0,x' - x) , 



where the fermion wave functions if)p r (t, x) are given by 



<r(*»x) 



(22) 



(23) 



with the Dirac spinors normalized to 

u A (p, r)u A (p, s) = 2m A 5 rs . (24) 

Using 3-momentum conservation at the vertices and, consequently, E^ A = E A , , the 
vertices in the last integral in Eq. (}22|) can be translated from times t < and t = to t = 
and — t > 0, just as we did in the scalar case [23] and also for the corrections to the vacuum 
energy above. Equation (|22jl can then be written in the covariant form 

(p A ,r\H BW \p' A , s) = (E v + E A A ) (27r) 3 5(p A - p^) 5 rs 

- ig 2 J d 4 xd 4 x'5(x°) [j$ Atr (x)S$(x - x')^r AjS (x')j A F (x - x') , (25) 

where the integral over d 3 x serves to implement 3-momentum conservation. 
Finally, then, the matrix elements of -f^ew take the form 

(p A , r|#Bw|PA, s) = (Ey + E$ A ) (2n) 3 5(p A - p' A ) 5 rs 

1 



+ 



2E L 



u A (p A , r)G(p A )u A (p A , s)} (27r) 3 5(pA - p' A ) , (26) 



where G(p A ) is written in momentum space by use of the non-covariant and covariant ex- 
pressions (114)) and ()13|) for the Feynman propagators in momentum space in Eqs. (|22j) and 
()25|) . respectively, to give the following equivalent expressions j2H] 

G(p) = -9 



-9 



d 3 p' 


1 


EpHo - p' -l + m A 


(2tt) 3 


2E$ 2^ p _ 


-p' E^+Up^-E^ 


d 3 p' 


1 


-Ep,iQ - p' ■ 7 + m A 


(2vr) 3 


2E$ 2^ 


p/ E$ + Up_p, + E£ 


d 4 p' 




p' ■ 7 + m A 



(2tt) 4 [(p - p'f -fi 2 + ie] \p' 2 -m\ + ie] 



(27) 



Po=E£ 
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The equivalence between the two expression is shown directly in Appendix [X] by performing 
the integration over p' Q in the manifestly covariant form. 

We will now show that the contribution to -ffew m the second line of Eq. (J26j) can 
be absorbed into a renormalization of the mass m A . To this end, consider the covariant 
expression for G(jp) in Eq. (}2"T|) : from covariance arguments, one has that G(p) is of the 
form 

G(p) = [G^)p ■ 7 + G (p 2 )m A ] po=E , , (28) 
hence by use of the Dirac equation 

u A (p, r)G(p)u A (p, s) = u A (p, r) [G^rrfyrriA + G (m A )m A ] u A (p, s) 

= 2m 2 A [G 1 (m 2 A ) + G (m 2 A )]S rs . (29) 

Consequently, we define 

Am 2 A = 2m 2 A [G x (m A ) + G {m A )] (30) 

in order to write 

Ev + E$ A + ^ (2vt) 3 5( Pa - p^) S rs . (31) 

We are now in a position to perform the mass renormalization, completely within the 
Hamiltonian framework. First, define the renormalized or physical mass through 



iA 

PA ^ O TP A 



= E V + M A , (32) 

p A =0 



then 



/\ TY) ~ 

M A = m A + + 0{g 4 ) . (33) 

IVTiA 

We can use Eq. (JHHj) to express uia in terms of Ma in Eq. (|3ip. in particular in -Ep 4 , taking 
into account that Am 2 A is of order (? 2 . Working consequently to order g 2 , we arrive at 

2 



*V + E$ A + ^ = E V+X /M 2 A + P 2 A + Otf) . (34) 

Remarkably, through the mass renormalization (to the order presently considered) we have 
obtained an expression for the energy which is exactly covariant, in contradistinction to Eq. 

en. 

Now, the arguments to arrive at the manifestly covariant expression in Eq. ()27|) (see 
also Appendix El) and from there to the crucial equation (}2"8*|) . are formal in the sense that 
the (unregularized) expression for G(p) is UV divergent. We show in Appendix [HI that a 
careful derivation using different regularization schemes leads to the same results. Observe, 
again, the absence of IR divergencies in the limit \x — >• in the expressions (|27|) (or in the 
corresponding explicit expressions presented in Appendix |B|). 
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2.3 Two-fermion states and effective potential 



Finally, in order to obtain the effective Schrodinger equation for AB bound states, we con- 
sider the Bloch- Wilson Hamiltonian for the subspace 



Q = span { \p A , r; p B , s) \p A , Pb G 



1.2} 



(35) 



of all (free) states of one fermion A and one fermion B, non-covariantly normalized as in Eq. 
(121)1 . Then the matrix elements of the Bloch- Wilson Hamiltonian to lowest non-trivial order 
turn out to be 

(Pa, r; Pb, s\H BW \p' A , r'; p' B , s') 

= (e v + ^M\ + p\ + yjMZ + plj (27r) 3 5(p A - p' A )6 rr ,(2n) 3 6(p B - p' B )S ssl 



dte~ eltl / d 3 xd 3 x' 







[ dte- eltl 


1 d 3 x d 3 x' 


/-oo J 





x A F (0-t,x-x') j$ A7r (t,x!)il$ Jt,x! 



^, r (0,x)^ r ,(0,x) 



x A F (0-t,x-x') ff B . a (t,x!)rf B s(t,x! 



(36) 



where Ma and M B are the renormalized masses defined as in Eq. (|32)h and the wave functions 
ipp^(t, x) have been introduced in Eq. (j23J). 

Due to the non-conservation of the perturbative energies, generally E^ A + E^ B ^ E^ + 
Ep, for pa + Pb = Pa + p'b, an d the corresponding lack of time translation invariance, the 
expressions with Feynman propagators in Eq. (jSSj) cannot be converted to on-shell Feynman 
diagrams. It has not been much emphasized in that the mixing of states with different 
perturbative energies is essential for the formation of bound states, because it is imperative 
for the localization of the constituents in relative position space: consider, e.g., states with 
total momentum Pa + Pb = 0, then a continuous superposition of states with different 
relative momenta p = p^ is necessary to obtain a wavefunction of finite extension in relative 
position space. 

Equation (j2HJ) leads to the following effective Schrodinger equation, 



Ma + Pa + \ M 2 B + p| <f>(p A , r; p B} s) 



d p A d p B . . J , , , i\ i/i ii i 

(Pa, r; ps, s\V\p A , r;p B ,s }<p{p A , r;p B ,s 



l%,J (2vr) 3 (2vr) 



(E-E v )cj>(pA,r;p B ,s) (37) 
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for the two-particle wave function in momentum space 

<P(pA,r;p B ,s) = (p A ,r;p B ,s\(j)) , \(j))eQ . (38) 

Through the use of the non-covariant representations ()14j) of the Feynman propagators in 
momentum space, the effective potential can be written as 

(PA,r;p B ,s\V\p' A ,r';p' B ,s') 

9 2 1 / I 



2EB 2E$ 2E% 2 ^pa-p' a \ E$ A + uj pa ^ a - E$ E$ B + u Pb ^ b - E* 

x [u A (p A , r) u A (p' A , r ')} [ub(pb, s) u B (p' B , s')] (2n) i 5(p A + p B - p' A - p' B ) . (39) 

The masses m A , B appearing in the kinetic energies E p ' b in the potential term can be replaced 
by their renormalized counterparts Ma,b to the present order in the perturbative expansion. 
Eq. (|39|) differs from its scalar analogue only through the products of Dirac spinors (cf. Ref. 

m)- 

The effective Hamiltonian commutes with the total 3- momentum operator P = + p B , 
hence we consider total momentum eigenstates from now on. In particular, we specialize to 
the center-of-mass system P = where the effective Schrodinger equation becomes 



2 

Um\ + + ^M B + pA 0(p; r, s)-g* £ f 

V ' r's'=l J 



d 3 P ' 



( 27r ) 3 J2E£2EZ2E*2E* 



+ 



2u p _ p , yE£ + u p _ p , - E$ E$ + ^ p _ p , - El 

x [u A (p, r) u A (p', r')] [u B (-p, s) u B (-p', a')] <j)(p'; r', s') — (E — E v ) 0(p; r, s) , (40) 

with the relative momentum p = = — p B and 

0(pa, r; p B , s) = <fi(p A ; r, s)(2ir) 3 5(p A + Pb) • (41) 

Observe in the potential term in Eqs. (|39|) and (|4(J|) the square roots of the kinetic energies 
which are characteristic of the non-locality of the interaction, and the differences of energies 
in the denominators which are due to the retardation of the interaction. The latter lead to 
a further non-locality in the effective potential, and also introduce non-hermiticity in the 
effective Hamiltonian. 

For the discussion of the non-relativistic and the one-body limits in the next section, and 
also for the numerical solution of the effective Schrodinger equation, it is convenient to cast 
Eq. PU]) into 2-spinorial form. To this end, the Dirac spinors are expressed in terms of Pauli 
spinors, most conveniently in the Dirac-Pauli representation, 

xa(p,0 

ua(p, r) = \jE^ + M A | P « } . (42) 

EjTW A XA ^ r) 
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Here the Pauli spinors are normalized in the usual way, 

Xa(p, r )XA(p,s) = S rs . 



(43) 



The Pauli spinors may or may not depend on the momentum p. The possibility of a momen- 
tum dependence is important if one wishes to employ helicity eigenspinors. In the present 
work, however, we will not make use of a momentum-dependent basis. 
Eq. ()40j) can now be rewritten in spinorial form as 



Mj + p2 + WM| + p2U(p 



d 3 p' IE£ + M A EZ + M b E$ + M A E* + M B 



(27T) 



2E$ 



2E* 



2E$ 



2E* 



1 



x 



+ 



1 



2u p _ p , \ E£ + - Ei E* + cj d _ d , - E* 



'p-p 



p p 



p ■ a A p' <t a 



E£ + M A E$ + M A 



p • a B p • cr B 



E* + M B E* + M B 



{p') = (E-Ev)(f>{p), (44) 



where the spinorial wave function 0(p) is defined as 

0(p) = ^( p] r ' s "> \Xa{y>, r) ® Xs(-P, s)] 



(45) 



As usual, cta is understood to act on XA(p, r ) only. 



3 Limiting cases, identical fermions and antifermions 

In the present section, we will perform a series of non-trivial analytical checks on the effective 
Schrodinger equation. We will begin with the non-relativistic limit: if the wave function 
0(p; r, s) is strongly suppressed for p 2 > M 2 , where 



M, 



M A M B 
Ma + Mb 



(46) 



is the (renormalized) reduced mass, we can approximate the effective Schrodinger equation 
in the center-of-mass frame, Eq. by 



P 



/ fX 2^ ,^ (P') =(e-E v -M a - M B )m 
J (2ix) A /j, 2 + (p - p') 2 \ J 



2M r ryir/ J (2tt) 3 /i 2 + (j 
where we have introduced the "fine structure constant" 



a 



9_ 

47T 
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(47) 



(4£ 



(for details on the approximation of the energy denominators, see Ref. |23j). Observe in 
particular that the effective Hamiltonian acts trivially on the spin degrees of freedom in this 
limit, as expected from non-relativistic scattering processes where spin orientations remain 
unchanged. As a consequence, Eq. (|4*7j) has the same form as in the case of scalar constituents 

The non-relativistic limit of Eq. (|4"4*j) is hence precisely what we expect on physical 
grounds, a non-relativistic Schrodinger equation with the usual Yukawa potential (after 
Fourier transforming to position space). As discussed in detail in Ref. [23], the limit is 
attained for a <C 1 and /i <C M r , both conditions being necessary. 

We will now consider the so-called one-body limit "Mb — > oo" where M\ 3> M\ and the 
wave function 0(p; r, s) is negligibly small for p 2 > M\. In this case, we can approximate 
the effective Schrodinger equation (|4"4*|) by 



'M 2 A + p 2 



IP) " 9 



d 3 p' 
(27p 



<E£ + M A E$ + M A 



2EA 



2E$ 



2u. 



p-p' 



P-P' 



P • tr A 



P o-a 



E£ + M A E$ + M A 



P 



= {E-E v -M B )<f>{ip). (49) 

The effective Hamiltonian acts trivially on the spin of particle B in this limit. Eq. (|49J) has 
the form of a relativistic equation for fermion A in an external potential, independent of the 
mass (except for a constant shift in the energy) and the spin orientation of particle B, in 
accord with physical expectations [TU| . 

We can, however, go one step further and compare Eq. (|49J) with the relativistic equation 
for particle A in the external potential due to a fixed source which exchanges spinless bosons 
of mass /i with particle A. The latter physical situation can also be described within the 
same general formalism, providing an internal consistency check for the application of the 
generalized Gell-Mann-Low theorem. 

To this end, we begin by defining the Hamiltonian H' of the fixed source system, which 
consists of the Hamiltonians corresponding to free fermions A and scalar bosons of mass /a, 
and the interaction term 



H[ = g d 3 x: ^(x)^(x)^(x) : +3^(0) • 



(50) 



The form of the second contribution to the interaction Hamiltonian results from replacing 
the dynamical fermion field B in Eq. (fTU|) with a fixed source at x = 0. To see that, we 
consider ^s(x) as a classical Dirac field with probability density given by 



P(x) 



4' 



x)^b(x) w ^ b (x)^b(x) 



(51) 



the latter approximate equality holding when ^>b( x ) describes a particle (not an anti-particle) 
and the relevant momenta satisfy p 2 <C M s . For a fermion B localized at x = we then 
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have 



Vb(x)Vb(x) = 5(x) , (52) 

from which Eq. (J50|) follows. 

The Bloch- Wilson Hamiltonian in the one-fermion sector is represented diagrammatically 
to lowest non-trivial order [see Eq. pi)) ] in Ref. (adding spin labels to the external lines 
for the present case). The corresponding algebraic expressions lead to a mass renormalization 
for fermion A which is identical to the one discussed in the previous section, and the following 
matrix elements of the effective Hamiltonian: 



*,x)V# Ait ,(t,x) 

A F (0-t,x-O), (53) 

where the (renormalized) vacuum energy E' v differs from Ey as determined from Eq. p2jl 
for the case of two dynamical fermions, precisely because fermion B has now been replaced 
by a fixed source which is formally considered to be part of the vacuum. 

Using the non-covariant expression (JT4"jl for the propagators and Eq. ()2*3*j) for the fermionic 
wave functions, Eq. leads to the effective Schrodinger equation 



(p A ,r\H BW \p A ,r) 

E' v + ^M| + p^ (2n) 3 5(p A - p' A )5 r 



-oo 




dte~ eltl / AA F (0-t,O-x) \"K A ,A 



-ig 2 / dte-^ / d 3 x ^ r (0,x)?/d r ,(0,x) 



Mi + p 2 0(p)-^y 



dV Ie£ + m a e$ + m a 



2E$ 



x 



2uj r ,_ r .i \ Ld, 



+ 



P P 



1 - 



p (T 



E£ + M A E* + M A 



P 



where the wave function is now defined as 

0(p) = ^2<P(P,r)x(P,r) = ^(p,r|0)x(p,r) 



{E-E' v )<j>{p), (54) 



(55) 



Equation (|54|) is identical to Eq. (|4^j) except for an irrelevant shift in the vacuum energy and 
the fact that the wave function in Eq. (|4"5j) includes the orientation of the spin of fermion 
B which, however, has no influence on the dynamics of particle A. Self-consistency of the 
method in the one-body limit is hence established. 
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We will now investigate how the effective Schrodinger equation changes when we replace 
one of the constituents, say fermion A, by the corresponding antiparticle. First of all, consider 
the one-A-antifermion sector where the matrix elements of the Bloch- Wilson Hamiltonian 
to lowest non-trivial order are given by 



[p A , r\H BW \p' A , s) = (e v + (2n) 3 5(p A - p' A ) S rs 



+ ig* 



+ ig 



CO 





dte~ eltl / d 3 xd 3 x' 



[compare with Eq. 
functions defined by 



^ „(t, *!)Sf(t - 0, x' - x)^ (0, x) A F (0 - t, x - x') 



dte~ eltl / d 3 xd 3 x' 



^, s (0, x)^(0 - t, x - x>pV (t, x') A F (t - 0, x' - x 



(56) 



and note the change in sign]. Here we need the antifermion wave 



<r(*> x ) 



2E£ 



(57) 



We can convert Eq. (j56|l to the form corresponding to a particle A by introducing the charge 
conjugate wave functions 



^f(t,x)=C tf£ r (t,x 



(5? 



(the superindex T stands for transposition), with the charge conjugation matrix C = —i'-f '-/ 2 
in the Dirac-Pauli representation. Eq. ()56|) is then transformed into expression for 



the subspace of one fermion A (including the sign) with the wave function ip^ r (t, x) being 



replaced by ippy(t,x), which corresponds to the replacement of the particle spinor wa(p, r) 
with the charge conjugate antiparticle spinor 



(59) 



for the description of the spin orientation of the antifermion. Since v A (p,r) is a positive- 
energy solution of the Dirac equation, the mass renormalization from m A to M A for an- 
tifermions is identical to the one for fermions A. 

We now turn to the two-particle sector with an antifermion A and a fermion B. The 
matrix elements of the corresponding Bloch- Wilson Hamiltonian to lowest non-trivial order 



17 



read 



(pa, r ; pb, s\h bw \p a , r '; Pb> s ') 

= (e v + ^M A + p A + ^M B + p B ^j (2tt) 3 5(Pa - p^)5 rr ,(27r) 3 5(p B - p' B )<W 



+ ig 2 



f 


dte~^ 


1 d 3 xd 3 x' 


1 — oo 






/* 


dte~ eltl 


1 d 3 x d 3 x 


/ — oo 







^ B ,,(0,x)^ >s ,(0,x) 



x A F (0 - t, x - x') x')^, r (t, x') 



^, ,(0,x)^ r (0 



x A F (0-t,x-x') ^ B>s (*,x')^y(*,x') 



(60) 



to be compared with Eq. (jHEJ). Again, the introduction of the charge conjugate wave func- 
tions converts Eq. (J5U|) to the form of Eq. (JSHJ), with the wave function i^p r (t, x) replaced 
by ipp'r(t, x), or the spinor -u^p, r) by f^(p, r). The effective Schrodinger equation for 
constituents Ai? can then be written in the spinorial form of Eq. (jjlj) , where in the defini- 
tion (|45)l of the wave function the Pauli spinor X^(P> r ) has to be replaced by the charge 
conjugate spinor 



(61) 



C,A{p,f) being the Pauli spinor that describes the spin orientation of the corresponding 
negative-energy solution of the Dirac equation. The most important result is that, just 
as for the mass renormalization, the interaction via a scalar boson is "charge conjugation 
blind", i.e., does not distinguish between a fermion A and its antifermion A. This is, in fact, 
the expected behaviour. 

The same arguments are used to describe an antifermion A interacting with a static source 
which leads to a Bloch- Wilson Hamiltonian and an effective Schrodinger equation analogous 
to Eqs. and (|54*jl. Consequently, the one-body limit is self-consistent also in the case of 
a (light) antifermion A. More interesting is the case of an antisource corresponding to the 
one-body limit Mb — > oo for an antifermion B: a simple-minded argument replaces ^b(x) 
in the interaction Hamiltonian with a classical negative-energy Dirac field, which leads to a 
probability density 



p(x) = ^Jj(x)V>b(x) « -V>b(x)V> 



B X 



(62) 



if we suppose that the relevant momenta satisfy p 2 <C M B . Equation (J52j) with p(x) = 8{x) 
would be in conflict with the one-body limit M B — > oo for constituents AB or AB. The 
reason is, of course, that p(x) is to be interpreted physically as a charge density (when multi- 
plied with the charge of fermion B), and is negative for antifermions, hence it is ■0b(x) , 0b(x) 
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which turns out to be positive and is to be replaced by <5(x). The latter results are related 
to the use of anticommutators in the quantization of the Dirac field. 

In order to have a formally satisfactory description, we define an approximately localized 
antifermion state with spin orientation s as 

|x = 0, Ax; s) = J |g| (|ax 2 ) ^ e-A-Va/s \p B , s) , (63) 

where |ps, s) denotes a 5-antifermion 3- momentum eigenstate. If we choose Mb large 
enough for 

Ax 2 M 2 B > 1 (64) 

to hold, we obtain 

/ 3 \ 3 / 2 

(x = 0, Ax; s\ : M*)M*) : |x = 0, Ax; s) = —— e -^ 2 /{^) (65) 

In the limit Ax — > the right-hand side of Eq. tends towards 5(x). Equation ()64|) implies 
that we need M B — > oo (even faster) in this limit. We take Eq. as justification to replace 
[: , 0b(x)'0b(x) :] with <5(x) for a fixed antisource, leading to the interaction Hamiltonian (}50|) . 
The one-body limit Mb — > oo is hence fully consistent also for the case of an antifermion B. 

We will now consider bound states of identical fermions, ^4-fermions to be concrete. To 
this end, we calculate the Bloch- Wilson Hamiltonian to lowest non-trivial order for the 
subspace 

Q = span{|p A ,r;p^,s)|p A ,p^ e M 3 ,r,s = 1,2} (66) 

of two-A-fermion states. As a consequence of the identity of the constituents, "crossed" 
diagrams appear which carry a relative minus sign due to the antisymmetry 

\PA,r;p' A ,s) = -\p' A ,s;p A ,r) . (67) 

However, mass renormalization works exactly as before, and also the effective Schrodinger 
equation is the same as in the case of AB bound states (with Mb = Ma) when we take into 
account the antisymmetry of the wave function 

0(Pa, r; p^, s) = -(j)(p' A ,s;p A ,r) . (68) 

In the center-of-mass system, we have consequently 

<Kp) = -<P(-PY (69) 
with the definitions (}4~Tj) and (|45|). where the transposition t refers to the tensor product, 

[xa(~P, r) ® xa(p, s)f = Xa(p, s) ® Xa(~P, r) . (70) 
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The consequences of antisymmetry for the solutions of the effective Schrodinger equation 
will be discussed in the next section. In the non-relativistic limit, however, the situation is 
particularly simple because the spin degrees of freedom do not participate in the dynamics: 
the solutions of Eq. (JTTj) with even orbital angular momentum have symmetric spatial wave 
functions and hence antisymmetric spin states (total spin zero), while solutions with odd 
orbital angular momentum have antisymmetric orbital wave functions and hence necessarily 
symmetric spin states (total spin one). 

In a theory which only contains A-fermions and scalar bosons, the results for the mass 
renormalization and the AA bound states are the same as the ones presented above for a 
theory with A- and 5-fermions, only the (irrelevant) corrections to the free vacuum energy, 
Ey — Eq, are different (the 5-fermion vacuum loops are missing). 

The last case we will consider in this section is the one of AA bound states, of one fermion 
and the corresponding antifermion. In this case, there are additional contributions from the 
virtual annihilation diagrams to the effective Hamiltonian, so that the effective potential 
reads 



(PA,r;PA, s\V\p' A ,r';p' A , s f ) 







I dte~ eltl 




' — oo 





V^, s ,(0,x)^ s (0,x 



A 

Pa. 







[ dte~ eltl 


1 d 3 xd 3 x' 


Loo J 





dte' €ltl / d 3 xd 3 x' 



x A F (0 - t,x - x') ^ ir (t,x')^(t,x') 



^ Air (0,x)^ ir ,(0,x) 

x A F (0-t,x-x') [^ s ,(*,x')< )S (*,x') 
^, r (0,x)^(0,x)_ 

x A F (0-t,x-x') jS ,(t,x')V^(t,x') 







[ dte~ eltl 


1 d 3 x d 3 x' 


Loo J 





V4y(0,x)Vd r ,(0,x) 



x A F (0-i,x-x') ^ 4 , r (t,x')^ iS (t,x') 



i 



+ 



i 



2E A A 2E£ A 2E$ a 2E A _ { ^pa-p' a \E£ a + ^ Pa - Pa - E A _ + u Pa . 



1 



2to- 



pa+pa 



^P A+ P A + Et + Ei 



x [u a (pa, r ) u A (p' A , r ')\ Mp'a> s ') va(Pa, s)] 
1 

+ 



Pa 



U P'a+P'a 



E A , — E A , 

Pa Pa 
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x [u A (p A , r) v A (p A , s)) [v A (p' A , s') u A (p' A , r')} \ (2tt) 3 5(p a + p A - p' A - p' A ) . (71) 



We will leave the discussion of the effect of virtual annihilation on the bound state energies 
in different theories and the related instability of these states for a future investigation. 

To close this section, we will briefly comment on the corresponding results in a purely 
bosonic theory, where the constituents are chosen to be charged scalar bosons. The consis- 
tency of the non-relativistic and one-body limits in this case has been shown in detail in 
Ref. [23] • If we substitute one or both of the constituents by antibosons, there are, com- 
pared to the fermionic case, no additional minus signs from anticommutation relations to 
take into account and, of course, no spinor structures, consequently the whole argument is 
much simpler than for fermionic constituents. The results are, however, finally the same: 
mass renormalization is identical for antiparticles and for particles, and the interaction due 
to scalar boson exchange is universally attractive and does not distinguish particles from 
antiparticles. 

As for a static antisource, consider the charge density (properly to be multiplied by the 
charge of boson B) 



for a classical negative-energy solution of the Klein-Gordon equation, in case that the relevant 
momenta fulfill p 2 -C M B . The probability density is the negative of p(x), hence we would 
replace 



in the interaction Hamiltonian for a localized antisource. A more formal argument proceeds 
in analogy with Eqs. for the fermionic case, with the result (fTH|) . The one-boson 

limit Mb —>■ oo is then consistent also in the case of an antiboson B. 

Finally, in the case of identical bosonic constituents, the effective Schrodinger equation is 
the same as for bosonic AB bound states, only that the wave function has to be symmetric 
under particle exchange in this case. Since there are no spin degrees of freedom in the 
scalar bosonic case, the spatial wave function has to be symmetric, hence only even angular 
momenta are allowed. 

4 Numerical solution 

In order to actually solve the effective Schrodinger equation in the form (j44j) . it is convenient 
first to separate off the angular and spin degrees of freedom. A direct numerical solution of 
Eq. (|4*4"|) would lead to numerical instabilities for equal and opposite momenta, due to the 
presence of a singularity in the integrand. The effective Hamiltonian is rotationally invariant, 
hence it is natural to consider total angular momentum eigenstates. To make contact to the 



p(x) = <f) B (x)i—(f) B (x) - (p B (x)i—0 B (x) 



2M B <f>* B (x)<f> B (x) 



(72) 




(73) 
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usual spectroscopy, we choose to couple first the individual spins to a total spin S and 
then couple this spin with the relative orbital angular momentum L to the total angular 
momentum J. The usual construction with Clebsch-Gordan coefficients yields simultaneous 
eigenstates of J 2 , J z , S 2 , and L 2 which we will denote as 25+1 3^m(p), P = p/IpI- Explicit 
expressions are given in Appendix IU1 

The Hamiltonian (}4"4"j) contains the helicity operators p ■ <x^ and p • erg. These operators 
are hermitian and unitary, and in particular 



The helicity operators are invariant under spatial rotations, however, they are odd under 
spatial parity transformations which maintain the spin directions unchanged. Since the 
Schrodinger equation ()44j) contains only even powers of helicity operators, the effective 
Hamiltonian is parity even (the intrinsic parities of the constituent fermions have no use 
in the present context, and we will not consider them in the following). 

We hence have the conservation of total angular momentum J and spatial parity (— 1)', 
but a priori not of relative orbital angular momentum I or total spin S. For given J, I can 
take the values J (for 5 = 0) and J, J± 1 (for S = 1), / = J, J — 1 being excluded for J = 0, 
S — 1. Taking into account the conservation of (—1)', we can hence conclude without any 
explicit calculation that the effective Hamiltonian may mix states with I — J, S — and 
with I = J, S = 1 on the one hand (we will call this "S-coupling"), and states I = J — 1, 
S = 1 and I — J + l, S — 1 ( "L-coupling" ) on the other. The effective Schrodinger equation 
will then decay into pairs of coupled one- dimensional equations. In the special case J = 0, 
neither of the two mixings is possible. For future use we remark that the use of helicity 
eigenstates is expected to diagonalize the effective Hamiltonian in the S-coupled sector, thus 
slightly simplifying the calculations, although there is no reason why L-coupling should not 
occur in this case. 

For the actual solution of the effective Schrodinger equation (|4*4"j) . we need explicit ex- 
pressions for the application of the helicity operators. Again, it is clear from the fact that the 
helicity operators preserve total angular momentum and change parity, that the application 
of a helicity operator maps S-coupled states to L-coupled states and vice versa. The explicit 
expressions, as well as a rather pedestrian way to derive them, are presented in Appendix 
mi Finally, the integration over the angles, i.e., over p', in Eq. (J44|) can be performed with 
the help of a partial wave decomposition combined with the spherical harmonics addition 
theorem, 




(74) 
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An 



V(p,p,cos6) = 



21 + 1 



ai(p, P ') Y im(i>)YL(f>') 



(75) 



1=0 



m=—l 



where p = |p|, 9 denotes the angle between p and p', and 



21 + 1 
2 



i 



dcos9 Pi (cos 9) V(p,p' , cos 9) . 



(76) 
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(Note that we have changed conventions relative to the ones employed in our earlier work 

The whole procedure outlined above can be carried through independently of the boson 
mass fi. In what follows, I will focus on the case of massless bosons \i = 0, as discussed 
in the introduction. The explicit form of the effective Schrodinger equation in this case is, 
again, given in Appendix O No additional difficulties are expected in the massive case in 
principle, even though the existence of a critical coupling constant changes the qualitative 
features of the spectrum. 

The (pairs of) one-dimensional integral equations can now be solved numerically. To 
this end, the equations were converted to (continuous) matrix form and the corresponding 
two-dimensional integrals approximated by finite sums over a discrete two-dimensional grid. 
The distribution of abscissas took the logarithmic singularity of the integrand and the long 
range of the wave functions in configuration space into account. We have approximated the 
solutions by a finite linear combination of an appropriately chosen set of basis functions, the 
same we had used before in the scalar case Two parameters that determine the shape 
of the basis functions were optimized variationally. The orthogonality of the basis functions 
could be retained numerically to 11 to 14 decimal places. Both energies and wave functions 
converged with increasing number of integration points and basis functions. However, in 
our experience [2H] convergence does not guarantee the correctness of the solutions if the 
choice of basis is inappropriate. For this reason, we have also checked the residual rj(p) of 
the solutions 4>i(p) defined as (E — ifew)0i(p) = The determination of the residual 

ri(p) and the "point- wise" convergence of the wave functions were limited essentially by the 
redundancy of the grid points (up to 400) with respect to the number of basis functions (up 
to 40). The analogue of the Gibbs phenomenon in Fourier series was, in the worst case, of 
the order of two percent. 

The numerical solutions (for fi = 0) are shown in Figs. ^ and [21 for two extreme mass 
ratios, Ma = Mb and Mb — > oo (with Ma fixed), as functions of the (Yukawa theory) fine 
structure constant a = g 2 /An. Between these two extremes, the eigenvalues for fixed a can 
be seen to vary smoothly with the mass ratio in Fig. |3]for a = 1. In all figures, the energy 
eigenvalues are represented normalized to twice the non-relativistic ionization energy in a 
Coulomb potential, M r a 2 , where M r denotes the reduced mass. 

In Figs. Hand0 it is seen that in the small-coupling limit a — > 0, the energies tend to 
the non-relativistic Coulomb values 

M r a 2 1 . . 

~— J, n=l,2,3,..., 77 

as expected for fi = 0. In particular, we observe the characteristic degeneracies in this limit. 
For instance, we expect and find that the following states tend to the same energy eigenvalue 
—M r a 2 /8 (principal quantum number n = 2) for a — > 0: the ground states in the L-coupled 
J = and J = 2 sectors (states 2 3 P and 2 3 P 2 in the usual spectroscopic notation n 2S+1 Lj) 
and in the S-coupled J = 1 sector (2^), as well as the first excited states in the S-coupled 
J = and J = 1 sectors (2 1 S , and 2 3 Pi) and in the L-coupled J = 1 sector (2 3 Si). To be 
precise, the lowest-energy eigenstates in the S-coupled J = 1 sector are linear combinations 
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Figure 1: The lowest energy eigenvalues (principal quantum numbers n 
the effective Hamiltonian as functions of the fine structure constant a 
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Figure 2: The lowest energy eigenvalues (n — 1 and n = 2) in the one-body limit M B — > oo. 
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Figure 3: The variation of the lowest energy eigenvalues with the mass ratio M r /M A between 
the extreme cases M r = M A /2 and M r = M A depicted in Figs. 1 and 2, respectively, for 
fixed a — 1. 
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of the 2 x Pi and 2 3 P 1 states, the coefficients depending on the mass ratio. In Fig. only 
three of these degenerate states are visible because there are always two curves lying on top 
of each other (see the discussion of the one-body limit below). 

Let us now discuss the equal-mass case of Fig. ^ in detail. The binding is weaker in 
this case than predicted by the non-relativistic formula. Compared to an electromagnetic 
interaction (exchange of photons) as in positronium, where we eliminate the contribution of 
the virtual annihilation diagram, the sign of the relativistic corrections is just the opposite 
(with the exception of the l 3 5*i state). Also, the ordering of the levels is different. 

In the case of equal masses, the effective Hamiltonian possesses an additional symmetry 
under the exchange of the fermions A and B, 



in terms of the spinorial wave function in the center-of-mass system (see Eq. (170(0 . This 
symmetry is implicit in our discussion of bound states of identical fermions in Section 3, 
where we saw that in the case of identical constituents the effective Schrodinger equation is 
the same as for AB bound states, hence it must possess exchange symmetry. The exchange 
parity for the angular momentum eigenstates is (— 1) . Together with the symmetry 
under spatial parity (—1)' (remember that we do not take the intrinsic parities of the fermions 
into account), this new symmetry forbids the S'-coupling, hence S becomes a good quantum 
number in this case. This can also be seen explicitly in the expressions for the matrix elements 
of the effective potential in the S-coupled sector, Eq. (|169|) in Appendix O for Ma = M B . In 
the case of identical fermionic constituents, the wave function must be antisymmetric with 
respect to particle exchange, hence, for J even, in the S-coupled sector only S = states 
are possible, while we have both I = J — 1 and I = J + 1 states in the L-coupled sector, 
in particular the coupling between these states remains. On the other hand, for J odd, we 
have only S — 1 states in the S-coupled sector, and no possible state in the L-coupled sector. 
Among the eight states shown in Fig. five are antisymmetric under particle exchange, 
namely the states l^So, 2 1 S'o, 2 3 P , 2 3 Pi, and 2 3 P 2 - We can see explicitly in these examples 
that the absence of S-coupling is necessary for the antisymmetry of the states. 

In the one-body limit Mb — > oo depicted in Fig. |21 the sign of the relativistic corrections 
changes for several states with respect to the equal-mass case. Within numerical accuracy, 
there are always two exactly degenerate states. The reason for the degeneracy can be seen 
from the effective Schrodinger equation in the one-body limit, Eq. ()49|): the effective Hamil- 
tonian is invariant under rotations which involve the spatial coordinates and the spin 
of particle A only, and independently under rotations of spin s# (which does not affect the 
dynamics). The total angular momentum of fermion A, = L + s^, is then a conserved 
quantity, and it is natural to consider simultaneous eigenstates of j\, ja,zi L 2 , and sb, z - 
Since ja — I ± 1/2 and spatial parity (—1)' is conserved as before, it follows that I is a good 
quantum number in this limit. 

As a further consequence, states that only differ in the value of ]a, z or sb, z are degenerate, 
and one may as well consider simultaneous eigenstates of J 2 , J z , ]\, and L 2 , where J = }a+ s b 
with eigenvalues J = ja ± 1/2. Now states which only differ in the eigenvalue of J 2 or J z 
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are degenerate. We can express the latter states in terms of the simultaneous eigenstates of 
J 2 , J z , L 2 , and S 2 that we are using in the numerical calculations. From the beforegoing 
discussion, we expect the absence of L-coupling and the degeneracy of one of the S-coupled 
states which is eigenstate of ]\ with eigenvalue ja — J + 1/2 and I = J, with an L-coupled 
state with J 2 -eigenvalue J + 1 and I = J (hence necessarily ja = J + 1/2). The other 
S-coupled state with I = J is an j^-eigenstate with Ja = J — 1/2, and is degenerate with an 
L-coupled state with J 2 -eigenvalue J — 1, I = J, and ]a = J — 1/2. Explicit expressions for 
the eigenstates in the different coupling schemes and their relations are given in Appendix 
O where we follow a different line of reasoning starting from the explicit expressions for the 
matrix elements of the effective potential in the S- and L-coupled sectors, Eqs. ()169|) and 

(HZ2D- 

These expectations are fully borne out in the results of the numerical calculations. Among 
the eight states calculated in Fig. |2J we expect and find the following to be degenerate: 1 ^So 
and l 3 «Si, 2 1 5'o and 2 3 S*i, 2 3 Po with one linear combination of 2 1 Pi and 2 3 Pi, and 2 3 P2 with 
the orthogonal linear combination of 2 1 Pi and 2 3 Pi. The coefficients of these linear combi- 
nations are predicted and found to be — a/1/3 and a/2/3, and a/2/3 and a/1/3, respectively 
(see Eq. (I18U|) in Appendix^, and correspond to eigenstates of with eigenvalues ja = 1/2 
and 3/2. Unlike for the Dirac equation with an electromagnetic Coulomb potential, states 
with the same ja but different /, here ja = 1/2 and I = (states 2 l So and 2 3 Si) or I = 1 
(states 2 3 Po and the first linear combination of 2 1 P\ and 2 3 Pi), are not degenerate. In 
addition, the ordering of the / = 1 states is opposite to the electromagnetic case. 

Figure El shows the smooth transition between the two extreme mass ratios Ma = Mb 
and Ma/ Mb = for fixed a — 1, and the appearance of the characteristic degeneracies in the 
limit Ma/ Mb — > 0. One would like to compare the coefficients of the mixing of the S-coupled 
states 2 1 P 1 and 2 3 P X against theoretical predictions. However, the diagonalization of the 
effective potential matrix (|169J) is generally not possible without solving the entire equation, 
i.e., it cannot be isolated from the p-dependence of the wave function. To the order a 4 of the 
first relativistic corrections, the diagonalization can still be performed analytically and turns 
out to factorize from the "radial" p-dependence. The results are presented in Appendix |TJ] 
and provide very satisfactory approximations to the corresponding results of the numerical 
calculations. In particular, the analytical results indicate that there is no hyperfine splitting 
to the order a 4 and Ma/ Mb, for the energy levels that are degenerate in the one-body limit. 
In fact, even for a = 1 no hyperfine splitting of the order Ma/ Mb is visible in Fig. |21 

5 Conclusions 

We have presented what appears to be the first consistent treatment of bound states in 
Yukawa theory. It is the result of a straightforward application of the generalized Gell-Mann- 
Low theorem. The consistency of the method has been checked thoroughly. In particular, we 
have shown that mass renormalization can be performed exactly as in a manifestly covariant 
formulation, even though the renormalization conditions were imposed entirely within our 
Hamiltonian framework. We have checked the non-relativistic and one-body limits, replaced 
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the fermionic constituents by antifermions, and considered the case of identical constituents. 
In all these cases, the formalism generates the correct results in a very natural way. In the 
numerical calculations, no abnormal solutions have been found (nor were there expected 
to be any, due to the absence of relative time or energy as a dynamical variable), and all 
the characteristic degeneracies in the non-relativistic and one-body limits show up in the 
numerical results with very good accuracy. 

In general terms, the framework presented here has several advantages over other formu- 
lations of quantum field theoretic bound state equations. As we have shown, the derivation of 
the effective Schrodinger equation is straightforward and presents no essential complications 
in the case of fermionic constituents as compared to scalar bosons. In principle, the complete 
bound state spectrum can be obtained as we have demonstrated by numerically determining 
the eight lowest-lying states (corresponding to the non-relativistic principal quantum num- 
bers n — 1 and n = 2). The wave functions for these states are also obtained in the course 
of the (approximate) diagonalization of the effective Hamiltonian. 

Several rather technical issues, which are nonetheless expected to be important for re- 
lated work in the near future, have been treated in detail in the appendices. In particular, 
the relation between manifestly covariant and non-covariant representations of the relevant 
loop integrals has been established (Appendix El . We have discussed dimensional, Pauli- 
Villars, Schwinger proper time and covariant and non-covariant cutoff regularizations for 
the appearing one-loop integrals from a non-covariant Hamiltonian perspective (Appendix 
IB]) . Finally, we have presented explicit expressions for the angular momentum eigenstates 
in different coupling schemes, for the application of helicity operators and the coefficient 
functions in a partial wave expansion, all necessary ingredients for the separation of angular 
and spin variables in the effective Schrodinger equation (Appendix EJ) . 

The results of this work, if only as a point of departure for the application to more realistic 
physical situations in the future, bear on fundamental issues in nuclear and high energy 
physics, as for example the nucleon-nucleon interaction. In this respect, one interesting 
particular result of the numerical computations is the qualitative difference between the 
relativistic bound state spectra for scalar (boson exchange) and electromagnetic interactions. 
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A Covariant and non-covariant representations of loop 
integrals 

In this appendix, we will relate different expressions for loop integrals in momentum space, 
where the manifestly covariant representations arise directly from the momentum space 
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Feynman rules, while the non-covariant representations result naturally from the application 
of the Gell-Mann-Low theorem. 

Let us begin with the lowest-order correction to the vacuum energy density Eq. (|19|) . 
The equivalence will be established by performing the integrations over po and p' through 
the use of complex integration theory. For greater clarity, we will first discuss the analogous 
problem in a purely scalar theory |23 r The corresponding formula differs from Eq. (j!9|) by 
the global sign (for both expressions) and the numerators of the integrands which are simply 
equal to one. Considering the first term (for particle A) for concreteness, we begin with the 
integral over p' Q , 



dPo 1 

2tx [(p — p') 2 — fi 2 + ie] [p' 2 — m 2 A + ie] 



(79) 



Standard application of the residue theorem, closing the integration contour through the 
usual large semicircle in the upper half plane, gives 
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This expression becomes much more transparent after a decomposition in partial fractions 
with respect to p , 
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where we have used the formula 



1 



UO — 17] OJ + it] 



27ii5(u>) 



^1) 



(82) 



(for 77 — > 0) . Equation (j81|) cannot be correct as it stands: the original integral (|79"Jl is 
even under p — > —p (by the substitution p' — * — p' of the integration variable), and this 
symmetry is manifestly broken by the delta function in Eq. ()81j1 . 



However, the term with the delta function only contributes for po 



Wp-p' - £y, and 

this is precisely the value of po where the two poles in the upper half plane coincide. The 
correct evaluation of the residue at the double pole gives for the integral in this case 



2u p _ p , 2E$ 



1 



Po + u p _ p/ + E§ - if] po- Wp-p' - E§ + it] 



— TP A 
P-P' p' 



so the result (j81|) is wrong for this value of po, and the term with the delta function has to 
be omitted in (j81j) . The integration over po can then be performed easily, using the residue 
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theorem again. The result, after substituting p' — > — p', is the one expected from Eq. (|T5j) or 
rather its scalar analogue. Note that the use of Eq. ()81j) including the delta function would 
lead to additional (incorrect) terms in Eq. (|19p. 

In the Yukawa case, Eq. (fTT?|) proper, the argument is nearly identical, only the ex- 
pressions are slightly more complicated. The result for the naive p -integration, after the 
decomposition in partial fractions, reads 



dp'p 

27T [(p 
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(84) 



By a calculation of the residue at the double pole for the special case p = u; p _ p / — E^,, one 
can again show that the term with the delta function in Eq. (J84|) is spurious. Integration 
over p of the rest gives the desired result. 

We now turn to the lowest-order corrections to the mass, Eq. (|27|). As before, we begin 
with the scalar case where all the numerators in Eq. (|27|) are replaced by one [2H]- Then the 
result of the pg-integration in the manifestly covariant expression is given precisely by Eq. 
()81|) above, with p to be replaced by E£. The delta function in Eq. (|5T|l is, of course, again 
spurious, although it cannot give any contribution anyway as long as /i 2 < 4m\. 

In Yukawa theory, where Eq. ()27|) properly applies, the pg-integration gives the following 
result after a decomposition in partial fractions (with respect to po), 
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putting po = -Ep . Again, the delta function turns out to be spurious (by explicitly considering 
the case of a double pole in the upper half plane), which leads to the non-covariant expression 
in Eq. (jZ3). 



B Regular izat ion of one-loop integrals 

The aim of this appendix is to derive Eq. (|28|) starting from the non-covariant expression 
for G(p) in Eq. (}2Tj) in a suitably regularized form, so that all integrals appearing in the 
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derivation are well-defined. By far the simplest way is to go through the manifestly covariant 
form also presented in Eq. ()27|). regularized correspondingly for the present purpose. For the 
complications that arise in a direct derivation in the non-covariant formulation for the simpler 
case of a purely scalar theory where it has to be shown that G(p) is actually independent of 
p (and only depends on p 2 \ po=E A = Tn A ), see Ref. |2*7j . 

B.l Dimensional regularization 

The technically simplest regularization scheme (although not the most natural one in the 
present context, see below) is dimensional regularization. The idea is to continuously change 
the dimension of space to smaller values where all the integrals are well-defined, to be save 
to spatial dimensions smaller than two, i.e., space-time dimensions D < 3. We can then 
establish the relation between the non-covariant and covariant expression in Eq. (J27|) for 
these dimensions and show that the dependence on D is analytical, with a simple pole 
appearing at D — 4. As a consequence, Eq. (|2*Hjl can be shown to hold true for space-time 
dimensions arbitrarily close to (but smaller than) four by analytical continuation. However, 
the definition of the integrals in arbitrary (non-integer) dimensions is made precise only in 
the (Euclidean) covariant formulation, which makes this form of regularization somewhat 
unnatural for the non-covariant expressions. 

In detail, we begin by establishing the relation between the non-covariant and covariant 
expressions for G(p) in Eq. (|27ji. but for (D — 1) spatial dimensions instead of three where, 
to begin with, D < 3 in order that all integrals are well-defined. The analogue of Eq. ([27|) 
in D dimensions can then be shown by integrating over p' in the covariant form exactly as 
detailed in Appendix El [see Eq. (|85|) and the following remarks]. The oP -1 ;/-integration is 
not touched in this process. With the covariant form at hand, we can introduce Feynman 
parameters in the usual way: 
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P0=E£ 



> f 1 j f dD( l xp-y + m A 

ig \ dx 



o 



(2tt) d [g 2 4- x{l — x)p 2 — x/j, 2 — (1 — x)m 2 A + ie l2 



Po=E£ 



(86) 



denoting the dimensionally regularized form of G(p) as G° R (p), where e = 4 — D. In the last 
step, we have shifted the integration variable to q = p' — xp and used that for the term q ■ 7 
emerging in the numerator, the integrand is odd. These manipulations are unproblematic 
as long as we stick to dimensions D < 3 where the integrals are well-defined. 
It is now convenient to decompose G^ R (p) in two parts in analogy with Eq. 



G? R ( P ) = [G™(p 2 )p ■ 7 + G^{p 2 )m A ] po=EA . (87) 
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However, all space-time vectors in this equation are still D-dimensional. Replacing p 2 
by m\ and Wick rotating to Euclidean space gives 
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The corresponding expression for is equal to Eq. (|55j) except for an additional 

factor of x in the integrand. It is these latter integrals which are rigorously defined for 
arbitrary, continuous values of D. They are actually defined in such a way as to make them 
analytical functions of D, with a simple pole at D — 4. We can use analytic continuation to 
define all the beforegoing integrals for 3 < D < 4, leaving the established relations intact, in 
particular Eq. ([87)1 . In the limit e — > or D — > 4, via the standard formulae of dimensional 
regularization, 
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plus terms which tend to zero in this limit. In Eq. (}89|) . k is the renormalization scale, and 
g is left dimensionless also for D^4. The corresponding expression for G^(m A ) is, again, 
equal to Eq. ()89j1 except for an additional factor of x in the integrand. 

After an integration by parts, the evaluation of the integrals over the Feynman parameter 
x is straightforward and yields 
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(for e — > 0). Observe the simple pole at e = or D = 4 and the absence of IR divergences 
for [A — > 0. 

The result for G^^m^) gives immediately the explicit dimensionally regularized expres- 
sion for Am 2 A in the purely scalar case considered in Ref. (where the coupling constant 
g has the dimension of mass). For Yukawa theory, we define Am A by 

Ami = 2m 2 A [G™(m 2 A ) + G^(m 2 A )} , (92) 
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so that 



u A (p,r)G» R (p)u A (p,s) = Am 2 A 5 r 



(93) 



from the analogue of Eq. (|29j) . where the Dirac equation in D dimensions is used. Explicitly, 
from Eqs. $Kj§ and (|9"Tjl . 
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(94) 



which can now be used to redefine m A in terms of M A as in Eq. (|33|) and the limit e — > be 
taken with the (finite) physical mass M A held fixed. 



B.2 Pauli-Villars regularization 

We now consider Pauli-Villars regularization which turns out to be the most natural regu- 
larization scheme in the context of Hamiltonian (non-covariant) perturbation theory. It is 
effected by subtracting from the expression in Eq. (|2"7J). denoted for the time being as G(p, /x) 
to make its dependence on the boson mass /i explicit, a similar contribution for a ficticious 
"heavy" boson of mass A to define 

C(p)=G( P ,/i)-G(p,A). (95) 

It is then easy to verify by power counting that G^ v (p) is UV finite, both from the non- 
covariant and the manifestly covariant expression. 

The pg-integration in the covariant expression can then be performed again as in Ap- 
pendix El only that the following integration over d 3 p' is now well-defined. This establishes 
the equivalence of the non-covariant and manifestly covariant expressions for G^ v (p) [i.e., the 
regularized form of Eq. 1)27)1] rigorously. Starting from the covariant form, we now introduce 
Feynman parameters and shift the integration variable as in Eq. (|86j) to obtain 



G{p^) = ^ g 2 f l dx r A xp-l + mA 
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and the analogous expression for G(p, A) [remembering that its difference G^ v (p) is well- 
defined in 4 dimensions]. From Eq. (J96j) and the corresponding expression for G(p, A), it is 
clear that G^ v (p) is of the form described in Eq. ()28)l. The mass renormalization can then 
already be effected at this stage as detailed following Eq. (J2EJ), reading G^ v (p) instead of 
G(p). 

To evaluate G^ v (p) analytically, it is easiest to apply dimensional regularization to 
G(p, fi) and G(p, A) separately and analytically continue the result for the difference G^ v (p) 
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back to four space-time dimensions. We start with the form obtained in Eq. (|96|) above, re- 
place p 2 \ po=E A by m A , and separate G^ v (p) in two parts G\ w A (m 2 A ) and Go^(m\), analogous 
to Eq. ()28j) . The continuation to D space-time dimensions and a Wick rotation to Euclidean 
space leads to the expression (jSHj) for the contribution Go(m 2 A , /i). The integrations over 
d D qE and the Feynman parameter then yield the results given in Eqs. (}9T)j) and (j9~Tj) for the 
contributions G (m A , (i) and Gi(m A , fi), respectively. The expressions for G (m A ,A) and 
Gi(m\, A) are obtained by replacing /x — > A. In the (relevant) limit of large A, one analyti- 
cally continues the last terms in Eqs. (|9T)|) and (f9~TJ) to A 2 > Am 2 A and expands the complete 
expression in powers of m 2 A / A 2 (up to second order). This gives, after several cancellations, 



G {m A , A) 
£1(777/4, A) 



r I 2 , , x , A 2 

' 1 - - 7b + ln(4vr) - In — + 1 



(4vr 



(47T) 



- — 1e + ln(47r) 



1 A 2 1 
2 ln ^ + - 



(97) 



for A 2 ^> m\. The same results can be obtained somewhat easier by replacing (1 — x) 2 m\ + 
xA 2 — > xA 2 in the argument of the logarithm in Eq. (}89|) [and analogously for Gi(m 2 A , A)] 
from the start. 

We hence obtain, finally, for Go(m A ) and Gi{m 2 A ) in Pauli-Villars regularization 
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/x 2 / (4m 
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(98) 
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(A 2 ^> m 2 A ). We obtain a logarithmic UV divergence with A — > 00 and, again, the absence of 
IR divergences for fi —>■ 0. We emphasize that we have used dimensional regularization only 
as a convenient calculational tool here, and that G^\{ m A) an< ^ G\\{m 2 A ) are well-defined in 
four space-time dimensions from the start. 

The result for G*q ^(m 2 ^) gives directly the expression for Am 2 A in Pauli-Villars regular- 
ization for the purely scalar case of Ref . j2S] • For Yukawa theory, we have from Eq. (J30j) 
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(l0l) 

B.3 Schwinger proper time regular izat ion 

The other two regularization schemes that we will discuss, Schwinger proper time and mo- 
mentum cutoff regularization, are only effective in Euclidean space. Hence we have to write 
the non-covariant expression in Eq. (|77jl in Euclidean space in order to regularize. After 
evaluation, it can then be analytically continued to physical values of the external variables. 

Let us begin with the diagrammatic expressions in Eq. ()22|) . We rewrite the first con- 
tribution to the mass renormalization by use of the covariant representations ()13|) of the 
propagators and Eq. (}23|) for the fermionic wave functions as 



dte- €ltl / d 3 xd 3 x' 



V4, r (0, x)S#(0 - t, x - x')^, s (t, x') A F (0 - t, x - x') 



d 3 p' d 3 k 
(2vr) 3 (2tt) 3 



d 3 xd 3 x' 6 iPA X f dp '° MPA,r)(p' ■7 + m A )u A (p' A ,s)]e ip ' <x x ' } 



2E L 



2tt 



Pq 2 ~ P /2 — m A + ie 



X 



e ik-(x-x') 



e ip A - 



dko 

2tt k* - k 2 - + xe ^ / 2jB a 



dt e et e i ( p ' 0+ko ~ Po ' >t 



(102) 



P0=M 



In the latter expression we Wick rotate in the mathematically positive sense 



IE 



Po -»• Po = l Po , 



kn 



(103) 



where p' E and fcjf are real after the rotation. The sense of the rotation is determined by the 
position of the poles in the integrand (by the ze-prescription). Since p' Q E and k E take both 
positive and negative values, it is imperative to Wick rotate t, too, keeping the exponents 
imaginary in order that the integrand do not blow up. This implies the rotation 



t->t= -it 1 



(104) 



in the negative sense, where — oo < t E < 0. Keeping p fixed at p = Eg, would then lead 
to a divergence in the t^-integration, so that we have to rotate, in addition, 



Po -> Po = Wo 



(105) 



After performing the integration, we then analytically continue the result to p E = —iE^, . 
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As a result of this Wick rotation, the expression (j!02|) is written as 



d 3 p' d 3 k 
2E A 2E A , J (2tt) 3 (2vr) 3 

PA p' A 



d 3 xd 3 x'e-' ipA - x 



dpf [u A ( PA , r) (pf^ + p' ■ 7 s + m A ) u A (p' A , a)] e^*"*') 



2tt 
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rll,E _ik-(x— x') rO 

rffc Q g , pVx . / E E t E i(p ,E +k E_ p E )t l 



(106) 



pE ^~ iE t' A 



where we have defined e £ = — ie to assure the convergence of the ^-integration. Alterna- 
tively, we can assume e to have a positive imaginary part from the beginning. As for the 7 
matrices, we have chosen the common convention where 



Polo + p' • 1 E = P' ■ 1 , 



(107) 



i.e., 7(f = i7o and 'jf = y E = —7*. Finally, for the momentum-space representation, 
we integrate over d 3 x, d 3 x' and d 3 k (taking advantage of the three-dimensional 5-function 
resulting from the x- and x'-integrations). The resulting contribution to G(p) [compare with 
Eq. (J2EJ)] is 
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2 f dy f dp' E pjf 7f + p' • -y E + m A 



(2tt) 3 J 2tt (p' E ) 2 + p' 2 + m\ 



dk E 



dt E e €EtE e i{p '° E+k °- p ° )tE 



2tt (k E ) 2 + (p - p') 2 + [i 2 
Turning now to the second contribution to the mass renormalization in Eq. 
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(109) 



po=e£ a 



and going through the steps following Eq. ([102)1 . we arrive at the following contribution to 
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G(p): 



d?p' f dp' E p' E ^ E + p' ■ j E + mi 



(2tt) 3 ./ 2tt ($f ) 2 + p' 2 + 



dfvQ 



2tt (fc^ ) 2 + (p - p') 2 + /i 2 



(110) 



It is easy to check the correctness of the results (J1U8)) and (jllUj) by performing the integrations 
over p' E and k E (using the residue theorem) and finally over t E , leading to the non-covariant 
expression in Eq. (|2*7I) . 

After this preparation, we are now in a position to introduce the Schwinger proper time 
regularization by replacing 

and analogously for the other (scalar) propagator. We emphasize that only in Euclidean 
space the modification of the lower limit for the parameter integration corresponds to the 
exponential suppression of the propagator for large momenta. Consequently, the proper time 
regularization of the first contribution Eq. (jl08j) reads 



-9 



d 3 p' 

(2tt) 3 Jl/A* 
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pZ^-iEA 



The integrations over p' E and k E in Eq. ()112|) are readily performed. After shifting the 
Euclidean time variable t E — > t' E = t E + 2ia(3p E / '(a + j3), one may analytically continue p E 



to — iE£ directly in the integrand to obtain 
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3„/ roo 
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A little algebra shows the first exponent in this expression to be negative. Integration over 
t' E gives 
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(114) 
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where the complementary error function is defined as 

erfc(x) = — / dte~ t2 = ^ dte~ t2 . (115) 

V 71 " A V 71 " ./-oo 

For the second contribution to the mass renormalization we start from Eq. piOjl . imple- 
ment the Schwinger proper time regularization as in Eq. (jlll|) and integrate over p' E , k E , 
and t E (after a shift and the continuation p E — > —iE^) to arrive at 
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(27T) 
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The results (I114|) and (|116J) represent the proper time regularization of the non-covariant 
expressions in Eq. (|27jl . When we add them up, we arrive at the proper time regularized 
equivalent to Eq. (|2*7j) . 
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where we have used that [see Eq. ()115|) ] 
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dt e"* 2 
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Next, we show that the sum (|117|) of the non-covariant expressions coincides with the 
Schwinger proper time regularized equivalent of the covariant form in Eq. ([27)1 . To this end, 
first perform a Wick rotation to rewrite the covariant expression in Eq. (J27)) in Euclidean 

space, 

d A p' E 



G(p) 
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p • 7 + m,4 



(2vr) 4 [(^ - p' E ) 2 + /i 2 ] [(p /i? ) 2 + m] 



(119) 



where all the products of Euclidean 4-vectors (including squares) are understood to be Eu- 
clidean scalar products, i.e., taken with the positive Euclidean metric. The proper time 
regularization is introduced as in Eq. (Jill)) and yields 



G7(P) = ~9 2 Q 2 dad ^S^ H E lo +P'-1 E + m A ] 



x e 
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• (120) 
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Integrating over p' E and continuing p E to —iE^ results in Eq. (jll7|) . thus establishing the 
equivalence of the non-covariant and covariant proper time regularized expressions. 

Note that there is a much faster way to establish this equivalence by adding up Eq. (J 112)) 
and its counterpart from Eq. ()110|) . substituting t E — > —t E in the latter expression, and 
using the formula 



/oo 
dt E e-^e l " tE =2n5(u). (121) 
-oo 



This procedure directly leads to the form ()12U|) . so it indeed represents a considerable short- 
cut. However, we were interested in writing down the proper time regularized equivalent of 
Eq. ()27|) . which is why we followed the calculation through to expressions ()114j) and ()116|) . 

To establish the form ()28|) . we start again from the covariant expression (|12(J|) . but inte- 
grate over all four components of p' E this time with the result 



(4vr) 2 J, 



g f (I )' , _ m . ); ^ _ ,'j (122) 
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-am 2 A -/3fj, 2 +al3p 2 /(a+l3) 



a + f3 

from where relation ()28|) can be established. Note that the exponent in Eq. (|122|) is manifestly 
negative for p 2 \ P0= Eg = m 2 A . 

In order to evaluate the functions Gq \{jn 2 A ) and G\\{m A ) in proper time regularization 
explicitly, we change variables (a,/?) to (p, a), where 

Mil P 

p = a + p , a 



a + 



0<a<l, p>P(<r,A 2 ) = nu«(^, (1 ■ (123) 

The integration over p can now be performed and leads to 

G P A T (P) = ~j£y 2 J o 1 da[ap. 1 + m A ] E 1 ([(1 - ofm\ + a^]P(a, A 2 )) , (124) 
where we have introduced the exponential integral function 



oo g— t roo (i~ xt 



E 1 (x)= J dt— = J dt—, x>0. (125) 

Using the expansion of E\(x) for small values of the argument we can approximate 
^ T (p)by 
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(4tt) : 



/ da[ap-j + m A } [ lE + \n ([(1 - a) 2 m 2 A + afi 2 }P(a, A 2 ))] , (126) 



all higher terms in the expansion being suppressed by powers of 1/A. For the sake of 
readability, the demonstration of this latter assertion is relegated to the end of this section. 
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The results ()9U|) and (|91|) and the integral of ln(A 2 P(a, A 2 )) which is elementary, can then 
be used to establish that 
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As before, Gp^m^) directly gives Am A in the purely scalar theory, while in the Yukawa 
case we have 
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Finally, we will give the demonstration of Eq. ()126j) which was relegated to this point. 
To begin with, let us note that the approximation ()126j) is not trivial because the function 
P(cr, A 2 ) takes large values for a close to the limits of integration, and 1. We further remark 
that the extension of the integration domain toO<cr< 1, p > 2/A 2 [compare with Eq. 
f!123|) ]. in which case the approximation corresponding to Eq. ()126|) is straightforward, leads 
to a different result. 

This being said, we start from Eq. (J 124)) and integrate by parts in order to get rid of the 
exponential integral function. The result can be written as 
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In the exponent in this expression, we keep the divergent terms (for a — > and a — *■ 1) and 
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expand the rest in powers of a and (1 — a), respectively, leading to 
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the higher terms in the expansion of the exponential of the finite terms being suppressed by 
powers of A, as we shall see shortly. First, note that we can, by way of the substitutions 



i and y = _L_ 

y 2a y 2(1-0- 



(132) 



in the first and second integral in Eq. (|131|) . respectively, express the contributions that 
contain the derivatives of In a and ln(l — a) in terms of the exponential integral functions 
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For small x, we have the expansions 



E x {x) = -7 £ - \nx + O(x) 
E 2 (x) = 1 + 0(xlnx) , 
E 3 (x) = l/2 + 0(x), 



(134) 



which lead to 
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in the limit of large A. In an analogous way, we can see that the higher orders in the expansion 
of the exponential of the finite terms (for a — > and a — > 1) in Eq. ()13()j) are suppressed 
by powers of A, taking into account that the expressions resulting from this expansion and 
the first terms in the integrals ()131|) are continuous bounded functions over the intervals in 
question, and that furthermore the higher orders in the expansion carry inverse powers of A. 
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As far as the first terms in the integrals ()131|) are concerned, the limit A — > oo can be 
taken naively there (note that the A-dependence in the argument of the logarithm is spurious 
since its derivative with respect to a gives zero), as we will show in a moment. If we suppose, 
for the time being, that this is indeed correct, the exponentials can be replaced by one, and 
an integration by parts yields 
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(136) 



for A — > oo (the A-dependence cancels between the two terms). Using Eq. (jl35j) and the 
results ()90|) and (|91|) for the remaining integration, we finally arrive at the expressions (|127j) 
and f)128|) . thus confirming the approximation ()126|) . 

Let us now show in detail that it is, indeed, correct to replace the exponential on the 
left-hand side of Eq. (I136|) by one. To this end, take the integrals on the left-hand side of 
Eq. (|136|) as they stand, perform the substitutions Eq. (|132|) and decompose the result in 
partial fractions to find, after a somewhat lengthy calculation, 
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where we have introduced the notations 
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For simplicity, we have written this expression for the case /x 2 > Am\ where the square roots 
are real. To compare with Eqs. (|127|) and ()128|) in the end, one has to analytically continue 
the result to smaller values of /i 2 . 
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The integrals in Eq. (jl37j) are well-defined with and without the exponentials. However, 
this is not true (in all cases) for the integrals over the partial fractions individually. It is 
clear then, that the exponentials can be considered as regulating factors for the integrals 
over the partial fractions. After summing up the individual results, the limit A — ► oo can 
safely be taken, thus removing the regulator. In the end, this is equivalent to replacing the 
exponentials by one in Eq. (jl36J) . The integrals in (j!37|) can also be explicitly calculated 
with the help of Eqs. ()133j) and (J 134)) in the limit A — > oo. As a result, Eq. (|136|) is recovered 
[using Eqs. (J9T?|) and (J9~T)) to evaluate the integral in Eq. (|136j) ]. 

These comments conclude the demonstration of the correctness of the approximation 
f)126j) and hence of the results (|127|) . (jl28J) and (jl29J) for the Schwinger proper time regular- 
ization. It is clear from the above that the proper time regularization will not be the method 
of choice in a Hamiltonian (not explicitly covariant) approach, considering the difficulties in 
establishing the regularized form of the non-covariant expressions and the evaluation of the 
corresponding (one-loop) integral as compared to the other regularizations discussed before. 



B.4 Momentum cutoff regularization 

At last, we discuss the probably simplest regularization scheme available, the use of a mo- 
mentum cutoff. As it turns out, it does not have very simple properties when used in the 
present context. We start with the non-covariant expressions written in Euclidean space 
after a Wick rotation in Eqs. ()108)) and (jllO)) . The momentum cutoff is implemented by 
restricting the integration domain (in Euclidean space) to (p' E ) 2 < A 2 , hence the complete 
regularized expression is 
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We can perform the integrations over k E (with the help of the residue theorem) and t E in 
Eq. ()139|) . separately for the two contributions, to obtain the result 
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Adding up the integrands in Eq. (|14Uj) gives 
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which obviously coincides with the cutoff regularization of the Euclidean version (jll9j) of 
the covariant form in Eq. (|27j) . thus establishing the equivalence of the non-covariant and 
covariant expressions in Eq. (J27|) in the momentum cutoff regularized form. One can use 
rotations in four- dimensional (Euclidean) space, taking into account the form of the integrand 
in Eq. ()141|) as well as the invariance of the integration measure and the integration domain, 
to show that G^ c (p) is in fact of the form of Eq. fl^g} . 

In order to obtain the three-dimensional form of the non-covariant expression in cutoff 
regularization, we start with Eq. (|14(Jj) and perform the p^-integration in 
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(with the d p'-integration restricted to p < A ). To this end, it is easiest to decompose the 
integrand in partial fractions with respect to p' E , for the first non-covariant contribution 




uJp-p>)lo + p' • 1 E + m A i^pTo + P' - 7 +m A 



Po 



pE 



Pf 



p-p- 



ho + P' • 1 E + m A ~ iE p'lo + P' • 7 + ™a 



E* + u 



Po 



P P 



(143) 



The decomposition in partial fractions for the second contribution in Eq. ()14(J|) can be ob- 
tained immediately from the one above by realizing that the two contributions are complex 
conjugate to each other (if we consider p E and the Euclidean 7-matrices as real, for example 
by using the Majorana representation). 

The Po^-integration is then straightforward, although the result is rather lengthy, indi- 
cating that the momentum cutoff would not usually be the regularization of choice in the 
present not manifestly covariant context, either. However, if we extended the p^-integration 
to the whole real axis while sticking somewhat arbitrarily to p' 2 < A 2 , we would obtain the 
following simple result for the sum of the two contributions: 
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n E - 
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(144) 



Eq. ()144|) results directly from the non-covariant expression in Eq. ()27|) by restricting the 
3- momentum integration to p' 2 < A 2 . It is hence, in a not manifestly covariant approach, a 
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very natural regularization scheme. Although it is apparent that four-dimensional Euclidean 
rotational invariance is broken in Eq. (|144|) and that hence G^ c (p) will not be of the form 
(128)) for any finite value of A, it is not clear a priori whether the form (|28|) could not be 
recovered in the limit A — ► oo. In the special case of a purely scalar theory, it was shown 
in Ref. [2Z1 that covariance is indeed reestablished for A — ► oo in the sense that G^ c (p) 
becomes independent of p in this limit (and only depends on the square of the four-vector, 

V% 0=e a = ™ 2 a)- 

In the following, we will consider the integration (J142)) of the integrand (J143)) and its 
complex conjugate counterpart arising from the second non-covariant contribution in the 
limit of large A and compare it with G^ c (p). To make this comparison mathematically 
precise, we divide the integration over the spatial momentum p' into the two regions p /2 < K 2 
and K 2 < p' 2 < A 2 with an intermediate scale K which fulfills p 2 , m% /i 2 -C K 2 <C A 2 . Then, 
in the limit A — > oo, we can approximate the integral over the first region, 



K 



d 3 P ' 



dpf 
2tt 



K 



d 3 p' 
(2^)3 



IE 



d^_ 

2tt 



(145) 



hence the integral over p' 2 < K 2 of Eq. (jl43|) and its complex conjugate, taking special care 
of the analytic continuation of the logarithms resulting from the p^-integration, tend to 



-9 



A 
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iE^ E + p ' . ^ + mA 



(2tt) 3 2E$ 2lu p _ p , E$ + 



p p' 
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d 3 p' 



1 iE^ + p' • j E + m A 



(2vr) 3 2E% 2cu p _ p , E$ + cu p _ p , 



ipo 



iEA 



(146) 



The first corrections to this result, of relative order (K / 'm a){K / A) [which can be suppressed 
in the limit of large A through a suitable choice of K, e.g., K = (m^A) 1 / 4 ], cancel among 
the two contributions. Incidentally, the second corrections, of order (K/rriA){K/A) 2 , also 
vanish. 



Now, for the other integration region K 2 < p 
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< A 2 , we have p 2 ,??^,// 2 -C p /2 ,A 2 . 



A 



lengthy calculation leads to the following result: the leading terms in a systematic expansion 
vanish for both contributions, the subleading terms cancel among the two contributions 
for the 7q -coefficients or as a result of the (p' — > — p')-symmetry of integrand, integration 
measure and domain for the 7 £ -coefiicient, and it is hence the subsub leading terms that give 
the dominant contributions which turn out to be 
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d 3 p' pf 7^ + p • 7 A ' + 2m A 
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+ 9 Z 



arcsin(|p'|/A) 



2|P'| 



|P 



'I - p /2 /A 2 
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d 3 p' 3pf7^ - p ■ j E yj - p /2 /A 2 
(2vr) 3 6^ A 



(147) 
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The first term in square brackets in Eq. ()147|) coincides with the dominant contribution 
from the part of G^ c (p) that stems from the second integration region K 2 < p' 2 < A 2 [cf. 
Eq. (|144j) ]. and hence, when added to Eq. (J146)) . combines to G^ c (p). The other terms in Eq. 
(1147)1 are readily integrated [2H] in the limit Kj A —* through the substitution y = |p'|/A. 
The final result is 



A d 4 p' E p' E ^ E + p' • 7*' + m A 



(2vr) 4 [(pf - pff + (p - p') 2 + fi 2 ][(p' E ) 2 + p' 2 + m\] 

+ j£y Q " ln 2 ) (i ^ 70 " p ' 7 + 2mA) " j£y T2 (3 ^ 70 + p ' 7) (148) 

in the limit A — > oo. 

The result ()148|) has two important consequences: first, together with the expressions 
(|157|) and (J 158)) for Gj£ (p) below, it provides an explicit expression for G^ c (p) in the limit 
of large A. Second, and more importantly, it shows that Lorentz invariance is broken, even in 
the limit A — > oo, through the last term in Eq. (j!48|) . For example, it can explicitly be shown 
that this term leads to a non-covariant contribution to the energies of the one-particle states 
by following the procedure detailed in Eqs. (|31 )) -([34 )) . It is then clear that the non-covariant 
cutoff regularization is not suitable for a fermionic theory. On the other hand, in a purely 
scalar theory, given by the m^-coefficient in Eq. (1148)1 , Lorentz invariance is recovered in the 
limit A — > oo, in agreement with the result in Ref. ^ZT\ which was obtained by a different 
method. 

Let us now return to the expression for G^ c (p) given in Eq. 1)141)1 . In order to calculate 
Go(m A ) and Gx(m A ) in cutoff regularization, one could use four-dimensional spherical coor- 
dinates for the <i 4 p /£; -integration (see also below), which would also explicitly confirm that 
G^f c (p) is of the form of Eq. ()28j) for any finite value of A. However, as long as one is only 
interested in the result for large enough A, it is much quicker in the present situation to in- 
troduce Feynman parameters and proceed in analogy with Eq. (jHU)) [in (D = 4)-dimensional 
Euclidean space and for the integration domain (p' E ) 2 < A 2 ]. We thus arrive at 







(27r) 4 {{ q E)2 + x (l _ x )(pEy + xfJ 2 + (! _ ]2 



(149) 



p 



with q E to be integrated over the four- dimensional Euclidean domain (q E + xp E ) 2 < A 2 (as a 
result of the shift of the integration variable to q E = p' E — xp E ). To evaluate this expression 
it is easiest to calculate the difference 

AGZ>(p) = Gf c (p)-G™ c ( P ) (150) 

between G A 4C (p) and 

-i M ja„e „„e „,E 



rM c M 2 f , rd*q E xp E ■ 1 E + m A 

Gu (p) = — q I dx I - — — „ 

J J \(q E ) 2 + x(l - x)( V E ) 2 + X{1 2 + (I - x)m 2 A f 



(151) 



P^-iEA 
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with the integration domain restricted to (q E ) 2 < A 2 [observe that in this case the first term 
in the numerator in Eq. (|149|) does not contribute because of the symmetry of the integration 
domain]. The G? 4 g £ -integral in G^, c (p) can be evaluated by standard methods. 

Incidentally, it would in principle be possible to define cutoff regularization via Eq. (j!51|) 
[instead of Eq. ([149)1 ] . a choice that would obviously fulfill Eq. fl2~%|) . too. The non-covariant 
expression corresponding to Eq. (jl39|) in this alternative regularization is 



^ C (P) 



dx 



A' J4„IE 
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dktQ 



Polo + P' ' 7 + m A 



(k§? x + ul_ v , x + (p' E f (1 - x) + (E$) 2 (1 
o 



dt E e eEtE e i( - p 'o E+k o-Po)t E + 



x 



dt E e eEtE e- i{p '° E+h °- pE)tE 



, (152) 



where the (i 4 p /£; -integration is now over (p' E — xp E ) 2 < A 2 . When comparing this expression 
with Eq. ()139|) . it becomes clear that this regularization is not too natural in the present 
non-covariant context. Still, we can integrate over k E (again with the help of the residue 
theorem) and tE in Eq. ()152j) . The result for the first contribution is 



q 2 I dx 
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A ' d A p' E 
(2vr) 4 
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1 . iP E o ~ p'o E )Vx~ {P E - P'o E ) 2 x + 3^_ p , x + 3 {p' E f (1 - x) + 3 (£#) a (1 



x 



^ P - P 'X + (p' E f (I - x) + (E$Y(1 



3/2 



(153) 



to be analytically continued to p E — > —iE^. The result for the second contribution can again 
be obtained by complex conjugation from the above (considering p E and the Euclidean 7- 
matrices as real). The integrations over both x and p' E in Eq. ()153|) look forbidding. However, 
the sum of the two non-covariant contributions is readily seen to give G^, c (p) as defined in 
Eq. (j!51|) . after shifting the four-momentum integration variable to q E = p' E — xp E . 

Returning to the explicit calculation of G^ c (p), we note that the difference AG^,(p) 
is an integral over the difference of the four-balls (q E + xp E ) 2 < A 2 and (q E ) 2 < A 2 . We 
introduce four-dimensional (Euclidean) spherical coordinates with the fourth axis oriented 
in direction of p E and the corresponding polar angle denoted as x- Then for the four-ball 
(q E + xp E ) 2 < A 2 the integration along the radius runs up to -R(x); 



R( x ) = A- x\p E \ cosx + 0(1/ A) 



(154) 



in the limit of large A. Furthermore, the denominator of the integrand in the region between 
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the two four-balls can be approximated by A in this limit. We then have 

AG^(p) = -, 2 fdxto sin 2 X f RiX \ r ^^Xlf + x(p E ^)+ mA ^ 



Jo A "A (2tt)*A* 

while the other components of q E in the numerator do not contribute for symmetry reasons 
(integration over the other angular coordinates). The integrations in Eq. (j!55|) yield, in the 
limit of large A, 

AG A)A ,(p) - j^y 2 ^~ - > ( 156 ) 

where we have used Ip^^f = p E ■ j E and replaced p E — > —iE^ = —ip°. In particular, for a 
purely scalar theory, the difference AG^^,(p) tends to zero for A — > oo. 

Now from Eqs. (jl5U|) and (|151|1 . the above result (jl56j) and the result of the integration 
over the Feynman parameter x in Eqs. (190(1 and (|91(1 . we finally get the explicit expressions 



~ 2A At 1 - T^T arcctg,/ ^ „^ , (157) 




(158) 

(A 2 ^> m 2 ^). Again, the result for G^(m A ) gives directly the expression for Am A in 
momentum cutoff regularization for the purely scalar case, while for Yukawa theory we have 



Ami = J 3 ln A^ + 7 _ 1 ^ _ /3 _ ^ l_ 

2m 2 A (47r) 2 |2 m\ 4 2m 2 A \2m A Am 4 A J m\ 




(160) 



C Separation of angular variables and spin 

In terms of the well-known eigenstates Xs,m s of total spin S = + s# (S = or S = 1), 
one has the following explicit expressions for the eigenstates of J 2 and J z with eigenvalues 
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J(J + 1) and M: 

%r(p) =Y JM (p) Xoo 

3 y J j-i,M(i>) 



1 



^2J(2J - 1) 



y/(J — M — 1)(J - M) Yj-i M+1 (p) Xi,-i 



3 yj M (p) 
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1 

" v/2(J + l)(2J + 3) 

- v/2(J - M + 1)( J + M + 1) Y J+ltM (p) xio 
+ a/(J - M + 1)(J - M + 2) y J+liM -i(p) xu 



(•/>!), 

V( J + M + 1)( J + M + 2) Y J+1M+l (p) X i,-i 



(161) 



To determine the action of the helicity operators on these eigenfunctions, the explicit 
spinor representation 



p a 



cos$ sin^e t<p 



(162) 



can be used, together with the following special instances of the spherical harmonics addition 
relation [SO] : 



cos$Y^ m (p) 
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/(/ - m + I) (I - m + 2) 
+ V (2/+1K2/ + 3) (163) 

It is then a straightforward, though somewhat lengthy exercise to derive the following 
formulas for the application of the helicity operators to the total angular momentum eigen- 
states: 



(p • va) l y J JM (P) = \j 2j^i M-i,m(p) - y %i,m(p) 



(p • »a) M-i,m(p) = %/(p) - %r(p) 



p • °a) %,(p) = -V^TT ^-m*(p) " V 2JTT 337jJ+1 ' A/(p 



(P • M + i,m(p) = -\/2JTT ^(P) - V 2JTI % ' (p) ' (164) 



and 



(P ' <tb) ^m(p) - -y 27^1 3 ^-i,m(p) + y %!,m(p) , 



(p • *b) : 'yj ,.uiP) = -\/2jTi 1;y ^ (p) " V 27TT 3 ^ m(p) ' 



(p • * B ) 3 y J JM (p) = -yj . uii>) - y 2TT1 3 ^j+i,m(p) 



(p • o-b) %!,m(p) = y %/(p) - y %*(p) • (165) 

In the special case J = 0, the states 3 3^j_i a/(p) an d ^ja/ (p) do n °t exist, and on the right- 
hand sides for the application of one of the helicity operators to l yj M {p) and 3 yj +1 a/(p)j 
only one term remains. 

By use of Eq. (J7H|) , where now 

2/ + 1 f 1 , „ „ , „ N 1 



2/ + 1 f 

a i{p,P) = — 2 — y dcos6Pi(cos 



2|P-P'| 
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(for n = 0), the effective Schrodinger equation Eq. decouples into pairs of coupled one- 
dimensional integral equations. For the S-coupled states, we introduce the wave function 



0(P) = S 0o J (pWm(p) + S #(pWm(p) • 
The effective Schrodinger equation for the coefficient functions becomes 
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(169) 
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On the other hand, with the wave function 

0( P ) = L <t>U(p) 3 y J j-i,M(p) + L ^M z y J J+ iAv) > 

the effective Schrodinger equation for the L-coupled states becomes 
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X 
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2J + 3 



(172) 



In the special case J = 0, s (f) , { (p) = L 0j_x(p) = 0, and the equations for 5 0q (p) and L 0j +1 (p) 
decouple. 

The functions a,i(p,p') are quite straightforwardly calculated from Eq. f)166|) . for any given 
value of I. The results for the functions which are relevant for the numerical calculations 
presented here (I < 3) are: 
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In the actual numerical calculations we did not rely on these explicit expressions, but rather 
used a quasi-algebraic method which increases speed and accuracy. To this end, the integrand 
in Eq. (j!66|) is written, after changing variables from cos 9 to x = |p— p'|, as a (finite) Laurent 
series in x around £y — Ep (or Ep 3 , — Ep 3 for the second term). This requires expressing 
P;(cos#) as a polynomial in x, a task left to the computer. Each term in the Laurent series 
has a known integral which is simply inserted. Any value of I can be handled easily by this 
method. 

In the one-body limit Mb — + oo, the matrix for the effective potential in the S-coupled 
sector, Eq. (jl69j) . tends to 
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J + l y/J(J+l) \ aj +1 (p,p') 

y/J(J+l) J J 2J + 3 

This matrix is diagonalized by the orthogonal linear combinations 

( S <PM \ = ( ~VJ \ ^J-1/2(P) 

V s rt(p) J V VTTT J V2TT1 

and 

( S <i>M \ = ( \ ^j+i/ 2 (p) 

V 3 <Pi(p) ) \ VJ ) V2J + T ' 

with the respective eigenvalues 
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aj(p,p') p p' aj_i(p,p') 
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+ 1^27+1 £^ + MaJ^ + M a 2J + 3 J' 
Alternatively, we can write the wave function as 

^(p) = S ^-i/ 2 b)j-i/ 2 ^(p) (178) 



.+i/2^m(p) = ^/ ^m(p) + y %/(p) , (180) 

anticipating the notation j A yf M (p)- l n the special case J = 0, of course, s <t>{{p) = 0, 
S <Pj_ 1 / 2 (p) = 0, and j-i/23^jm(p) does not exist, since there is only one state in the S- 
coupled J = sector. 

The matrix ()172|) of the effective potential for L-coupled states, on the other hand, 
becomes diagonal in the limit Mb — > oo, hence there is no L-coupling in this limit. The 
diagonal matrix elements tend to 
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Then, from Eqs. ffTTTHl) and (fT7T|) . 5 0j_ 1/2 (p) in Eq. (fTTHj) and L 0jIi/ 2 (p) = in E Q- 

(EOD, 

^(p) = Vtl /2 (p) ! Wii?(P). ( 182 ) 
fulfill exactly the same one-dimensional equation in the limit Mb — > oo (for J > 1), 

5 ^-i/ 2 (p) \ <? 2 r , , ,2 lE£ + M A E$ + M A 



A+p { l <p j j'-Up) ) W pp 



x 



aj(p,p') _ P P' aj-i(p,P') \ ( S 4> J j-i/ 2 (p') \ 

^27+1 E£ + M A E$ + M A 2J-1 J V L <f>jZ{ /2 (p') J 









-1/2C 








-1/2 ( 










.,( 









-1/2 (p) 
-i/ 2 b) 



Likewise, 5 0j +1/2 (p) in Eq. (jT7Rjl and L 0j+ 1/2 (jo) = in Eq. (fTTHjl . 



v+1/2 

both fulfill the equation 



2J + 1 + M A E$ + M A 2J + 3 H V25 /2 (p') 



x 



(183) 



P) = L ^/ 2 (P) 3 ^ 1 (P), (184) 



^8 )■ (i85) 



As a result, every state is (at least) twofold degenerate in the one-body limit. 

The states j-i/2D ; jm(p) an d j+i/2^jm(p) defined in Eq. IjlKUj). as well as j-i/23 ; j M 1 (p) = 
3 3^jm"(p) an d j+i/ 2 3^jm 1 (p) = ^jm^p); can be rewritten in terms of the eigenstates V^(p) 
of ji, jA,z, and L 2 (where j A = L + s A ), 



v i-i/2/ M 1 / - v / /-m + l/2F iim _ 1/2 (p) 

v m/2 r , 1 f ^ + m + l/2Vi/2(p) . 

(P) = /7TTTT , , s > ( 186 ) 

V2T+1 V v //_ m+ i/ 2 y /im+1/2 p 
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as 



J-i/2^m(p) = y ^j,Af+i/ 2 (P) ® ( x J + y—2j~ *Wi/ 2 (p) ® o 

v w/m / 7+M j_ 1/2 / \ / J-M w _ 1/2 /l 

j-i/2^ M (p) = y ^m+i /2 (p) ® ! J - y — ^— r /,M-i/ 2 (p) ® o 

. + i/ 2 ^m(p) = W 2(J + 1) ^m + i/ 2 (p) ® ! J - W 2(J+1) ^,M-1/ 2 (P) ® ^ o 



w + i n J-M + l vJ+1/2 /0\ / J + M+l J+1/a 

z+i/^i (p) = w 2(J+1) ywv»(p) ® ^ i J + y ^jttt 7 ' M - i/2(p) 1 o 



187) 



hence they are simultaneous eigenstates of J 2 , J z , and L 2 . 
Using the formulas 

(P-^ 1/2 (P) = -^(P)> 

(P'*)v£ 1 *(P) = -y£1£(P), (188) 

proven, e.g., with the help of Eqs. (j!62|) and ()163|) . and consistent with the fact that the 
application of helicity operators conserves angular momentum = L + and changes 
spatial parity (—1)', it is easily seen that for the wave functions (j!78|) and (j!82|) the effective 
Schrodinger equation in the one-body limit, Eq. (J49|) . is equivalent to Eq. (J183j) . while for 
the wave functions (fT79l and (fPMD . Eq. flU is equivalent to Eq. (fPSoT) . 

The coefficient functions ai(p,p f ) of the partial wave expansion in the one-body limit 
Mb — > oo are defined by 



az(p,p') = — / dcosePticose) o| _ _„ [ + - A ■ _„ ^ , (189) 



2Z+1 /- 1 

// a \< H I-',', ri H ) 

2|p-p'| VIp-p'I ^ + |p-Pl-^y, 

see Eq. ()49|) . The explicit expressions in the cases I — 0, 1, 2, which are the ones relevant for 
the numerical results presented here, are 

3 J(^-^)(p+p'-b-p'i) 



ai(p,p') 



4pp' j 2pp' 

| P 2 +P /2 ln P + P' P 2 +/-(^-^) 2 ]n / ^ + (p + j/)-^ 
2pp' |p — p'| 2pp' I i?^ 1 + \p — p'| — £y 
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02(P,P') 



5 J 3[2(p 2 + p' 2 ) - (J5£ - ^) 2 ](£ p A - Ef)(p + p'-\p- p'\) 



8pp' 



4p 2 p' 2 

VI31 QTo^2 i „/2> 



(E* - ^) [(p + p'f -\p- pf] 3[2(p 2 + p") - {E} - E$ 



A\2l 



Ap 2 p 



2^/2 



2pp' 



+ 



+ 



3(p 2 +p' 
4p 2 p' 2 



,/2n 2 



In 



p + p' 



|p — p'| 



3[p 2 +p' 2 

4p 2 p /2 



A\212 



hi 



'Ef + (p + pf)-E£ 
E£+\p-p'\-E$ 



(190) 



They coincide with the limit Mb — > oo of the corresponding expressions in Eq. ()173|) . 

The consistency of the one-body limit is hence fully established at the level of the effective 
Schrodinger equation after the separation of angular and spin variables. 



D Approximate diagonalization of the effective poten- 
tial matrices 



To order a 4 , i.e., for the lowest-order relativistic corrections, we can approximate 

p p p p 



Ep + M A 2M A 



E* + M B 2M 



B 



(191) 



and analogously for p' instead of p, in the potential terms of the effective Schrodinger equa- 
tions ()168|) and (j!71|) . Furthermore, terms containing 



P 



P 



P' 



P 



E£ + m M A E* + M B E$ + M A E* + M B 



(192) 



can be neglected. The matrix ()172j) in the L-coupled sector can then be approximated by 
the following diagonal matrix 



2J+ 1 

/ aj-i(P,P') 
2J - 1 

V 



L Vj J ^j^p') V/_ M+1 (p,p') 
L V/ +1 ,j- 1 (p,p') l Vj J +1 , j+1 (p,p') 





aj+i(p,p') 
2J + 3 / 



pp' pp' 



AMI 



4M 2 



/ aj(p,p') 



\ 



2J+1 








\ 



aj(p,p') 
2J+1 ) 



(193) 



In particular, to order a 4 there is no L-coupling. We hence expect, at least for moderate 
values of a, the mixing for the solutions in the L-coupled sector (with J > 1) to be rather 
small. 
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The matrix elements (jl69j) in the S-coupled sector can be approximated by 



S V J o(p,p') = aj(p,p')~ 
S V 1 J 1 (p,p') = a J (p,p')- 

s V J 1 (p,p') = s V 1 J (p,p') 



PP' PP' 
AM\ 4M| 

pp' pp' 
AM\ + 4M| 



j aj-i(p,p') | a J+i(P'P') 



2J- 1 



2J + 3 



(j + i) Q-z-ifapQ + j a J+ i(p,p') 



pp 



pp 



AM\ 4M b 



y/J(J+l) 



2J-1 2J + 3 

aj-i{p,p') a J+1 (p,p') 



2J - 1 



2J + 3 



. (194) 



This approximate matrix, somewhat fortunately, can be diagonalized through a p- and p'- 
independent transformation. After some algebraic labor, one obtains the following eigenvec- 
tors: 



/ 



; {+) 

J J,x 



V2 



1 



y/1 + AJ(J + l)x 2 



1 + AJ(J+ l)x 2 



1/2 



1/2 



\ 



V2 



1 + 



1 

/ l + 4J(J + l)x 2 



1/2 



1/2 



V 



(195) 



T + 4J(J+ l)x 2 J ) 
(J > 1); where we have expressed the mass dependence through the parameter 



x 



M 2 B -M\ 
M 2 B + M\ 



(196) 



The corresponding eigenvalues are 
2J+1 V ^ tP,P)- 2J+1 



pp 



AM 2 1 + 



x 



1 + 4J(J + l)x 2 \ aj_i(p,p') 



2J+ 1 



2J-1 



+ 



1 V / l + 4J(J + l)x 2 \ aj+ifapQ 



2J+ 1 



2J + 3 



2J + 



2J+ 1 4M 2 , 1 



.r 



1 - 



v/l + 4J(J + l)x 2 \ aj.x^p') 



2J+ 1 



2J- 1 
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' Vl + 4J(J+l)x 2 \ a J+1 (p,p') 
2J + 1 j 2J + 3 



(197) 



These formulas are valid for x > or M# > M^. Comparing these results with the approx- 
imate diagonal matrix elements for L-coupling, 

1 yj-y /\ _ aj(p,p') pp' 2 aj-i(p,p') 

2J-1 JJ {P,P) 2J+1 AM\l + x 2J-1 ' 

1 l v j + i ( f) _ aMp 1 ) _ p^_ J_ ^wfapQ ri98i 

2J + 3 ^ JJ 2J+1 4Mil + x 2J + 3 ' 

one notes the coincidence in the one-body limit x — > 1. For x close to one, Eq. (J197)) can be 
expanded in x around one. Consequently, the energy levels that are degenerate at x — 1 split 
for Ma M b by terms of the order of (at least) (Ma/M b ) 2 , which points to the fact that 
there is no hyperfine splitting of the levels in the strict sense, i.e., to orders a 4 and Ma/M b . 
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